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Abstract. The formalism allowing one to account for the effect of a finite space- 
time extent of particle production region is given. Its applications to the lifetime 
measurement of hadronic atoms produced by a high-energy beam in a thin target, as 
well as to the femtoscopy techniques widely used to measure space-time characteristics 
of the production processes, are discussed. Particularly, it is found that the neglect of 
the finite-size effect on the pionium lifetime measurement in the experiment DIRAC at 
CERN could lead to the lifetime overestimation comparable with the 10% statistical 
error. The theoretical systematic errors arising in the calculation of the finite-size effect 
due to the neglect of non-equal emission times in the pair center-of-mass system, the 
space-time coherence and the residual charge are shown to be negligible. 



PACS numbers: 03.65, 25.75, 36.10 



Finite-size effect on two-particle production 



2 



1. Introduction 

The determination, on a percent level accuracy, of the breakup probability of the 7r + 7r~ 
atoms produced by a high-energy beam in a thin target is of principle importance for a 
precise lifetime measurement of these atoms in the experiment DIRAC at CERN [1-4]. 
This experiment aims to measure the lifetime r 10 of the 7c + tt~ atoms in the ground state 
with 10% precision. As this lifetime of order 10 -15 s is determined by the probability of 
the annihilation 7r + 7r _ — > 7r°7r°: 1/t\q ~ |a[j — <2q| 2 , the DIRAC measurement enables to 
determine the absolute value of the difference af]— a\ of the s-wave isoscalar and isotensor 
Tin scattering lengths to 5%. This represents a factor of 4 improvement of the accuracy 
achieved in previous studies [5] , and is comparable with the precision of the most recent 
experiments E865 at BNL [6] and NA48/2 at CERN [7]. The former is based on a study 
of K e4: decays and yields the statistical error of 6% in a[j; this measurement essentially 
exploits other experimental data together with dispersion relations (Roy equations), 
the systematic and theoretical errors being estimated on the level of several percent. 
The latter studied the cusp effect at the tt + ti~ threshold in the distribution of the 2n° 
effective mass in —> n^^ir decays and yields |a° — a§ | with a few percent statistical 
precision and ~ 5% theoretical error. Both these measurements are in agreement with 
the preliminary DIRAC result based on ~ 40% of the collected statistics jl] as well as 
with the prediction of the standard chiral perturbation theory [8] . 

It should be stressed that the theoretical prediction for the difference a®— or] depends 
on the structure of the QCD vacuum. Thus, on the standard assumption of a strong 
quark condensate one has a[} — = 0.374±0.006 fm [8]. With the decreasing condensate 
this difference increases and can be up to 25 % larger [9]. The precise measurements of 
the 7T7T scattering lengths thus submit the understanding of chiral symmetry breaking 
of QCD to a crucial test. 

The method of the lifetime measurement is based on the production of 7r + 7r~ atoms 
in a thin target and subsequent detection of highly correlated 7r + 7r~ pairs leaving the 
target as a result of the breakup of a part of the tt + tt^ atoms which did not decay 
within the target [10J. Clearly, the breakup probability is a unique function of the 
target geometry and material, the Lorentz factor and the ground-state lifetime of the 
7T + 7r _ atom. The analysis shows that, to achieve the accuracy of 10% in the lifetime, 
the breakup probability, in more or less optimal conditions, should be measured to 4% 

There are two methods [2] - extrapolation and subtraction ones - which can be used 
to measure the breakup probability Pbr (or a combination of the breakup probabilities 
in different targets) defined as the ratio of the number iV]f of breakup atoms to the 
number Na of the atoms produced in the target: 

P br = N%/N A . (1) 

The extrapolation method requires the calculation of the number of produced 7r + 7r~ 
atoms Na based on the theory of the final-state interaction (FSI) in discrete and 
continuous spectrum [TOj QJ], [12]. This calculation, as well as the determination of 
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iVjf, is not required in the subtraction method which exploits the data taken on at least 
three different targets made out of the same material but consisting of a different number 
of layers of the same total thickness. However, this method needs a factor of 7 larger 
statistics [T] and cannot yield the required precision within the approved time-scale of 
the experiment DIRAC. 

The FSI effect on 7r + 7r~ production is sensitive to the space-time extent of the pion 
production region mainly through the distance r* between the tt + and tc~ production 
points in their center-of-mass (cm.) system. In [TO], only the Coulomb FSI was 
considered and the r*-dependence was treated in an approximate way, dividing the pion 
emitters into short-lived (SL) and long-lived (LL) ones. It was assumed that r* = 
for pion pairs arising solely from the SL emitters and characterized by the distances r* 
much smaller than the Bohr radius \a\ of the 7r + 7r~ system (a = —387.5 fm), otherwise 
r* = oo. 

The finite-size correction to such calculated number of non-atomic 7r + 7r~ pairs in the 
region of very small relative momenta in the pair cm. system, Q <C 1/r*, is determined 
by the three dimensionless combinations r*/a, fo/r* and fo/a of r*, a and the s-wave 
7r + 7r~ scattering length: / = ^(2dQ + al) rs 0.2 fm. Typically (r*) SL ~ 10 fm so that the 
correction is dominated by the strong interaction effect and can amount up to ~ 10%. 

Fortunately, due to a small binding energy e& ~ (m 7T a 2 )^ 1 , the finite-size correction 
to the production probability in discrete spectrum at r* <C \a\ is nearly the same as that 
in continuous spectrum at zero energy. Since the calculated number Na of produced 
atoms is approximately determined by the measured number of non-atomic 7r + 7r~ pairs 
and the ratio of weighted means of the finite-size correction factors corresponding to the 
production in discrete and continuous spectrum, the finite-size correction would cancel 
out in Na, up to 0((r*/a) 2 ) and 0(/o/a), provided we could measure the number of 
non-atomic pairs in the region of very small Q <C 1/r* [T3], [H]. 

At small values of Q and r*, the relative correction to the number of non-atomic 
pairs is positive (due to the effect of the strong FSI ~ 2f /r*) and changes sign at r* ~ 
10 fm (due to negative finite-size effect of the Coulomb FSI ~ 2r*/a). It appears that for 
r* < 20 fm the correction shows a quasi-linear behaviour in Q up to ~ 50 MeV/c, with 
almost a universal negative slope. For larger distances r*, the slope becomes positive 
and has a non-trivial Q- dependence. If the pions were produced at small distances r* of 
several fm, one could safely neglect the non-universal correction 0((r*/a) 2 ) and use the 
quasi-linear Q-dependence of the correction factor to interpolate to Q = 0. However, 
there is a non-negligible tail of the distances r* > 10 fm due to particle rescatterings and 
resonances (particularly, r* ~ 30 fm in the case when one of the two pions comes from 
the cu-meson decay). In the experiment DIRAC, the finite-size correction can lead to 
about percent underestimation of Na and - to several percent overestimation of N^- As 
a result, the corresponding overestimation of the extracted lifetime can be comparable 
with the 10% statistical error and should be taken into account. 

We discuss how to diminish the systematic error due to the finite-size effect on the 
lifetime measurement of hadronic atoms, using the correlation data on identical charged 
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pions (containing the information about the distances r* between the pion production 
points in the same experiment) together with the complete phase-space simulations 
within transport models. 

The formalism accounting for the finite space-time separation of particle emitters 
is also in the basis of the correlation measurements of the space-time characteristics 
of particle production - so called particle inter jerometry or correlation femtoscopy (see 
reviews [T5], [161 EH (HI EE l2Ql [21]). In fact, the femtoscopic correlations due to the 
Coulomb FSI between the emitted electron or positron and the residual nucleus in beta- 
decay are known for more than 70 years; the sensitivity of the differential decay rate to 
the nucleus charge and radius is taken into account in terms of the Fermi function which 
can be considered as an analogue of the correlation function in multiparticle production 
(see [22] for a discussion of the similarity and difference of femtoscopic correlations 
in beta-decay and multiparticle production). The femtoscopic correlations due to the 
quantum statistics (QS) of produced identical particles were observed almost 50 years 
ago as an enhanced production of pairs of identical pions with small opening angles 
(GGLP effect). The basics of the modern correlation femtoscopy were settled by Kopylov 
and Podgoretsky in early seventieth of the last century; they also pointed out a striking 
analogy between the femtoscopic momentum correlations of identical particles and the 
spectroscopic space-time correlations of photons (HBT effect), the latter allowing one 
to measure the spectral width of the light source as well as the angular radii of distant 
(stellar) objects (see [22] and references therein). Besides the space-time characteristics 
of particle production, the femtoscopic correlations yield also a valuable information on 
low-energy strong interaction between specific particles which can hardly be achieved 
by other means (see [20] and subsection 4.3.5). 

The paper is organized now as follows. In sections 2 and 4 we give the basics of 
the theory of two-particle correlations due to the FSI and QS effects. Particularly, the 
formalism and assumptions behind the correlation femtoscopy are discussed in section 4. 
Sections 5 and 6 deal with the one- and two-channel wave functions in the continuous and 
discrete spectrum. In sections 3 and 7, we apply the developed formalism to estimate 
the finite-size effect on the pionium lifetime measurement in the experiment DIRAC at 
CERN. The results are summarized in section 8. In Appendices A and B we consider 
the effect of non-equal times and derive the analytical expression for the normalization 
effect of the short-range interaction on the wave function of a hadronic atom, modifying 
the usual n~ 3//2 dependence of the pure Coulomb wave function on the main quantum 
number n. The reader interested mainly in practical application of the formalism to the 
lifetime measurement of hadronic atoms can start reading from sections 3 and 7 and 
consult the rest of the paper to clarify the eventual questions. 

2. Formalism 

The production of two particles at small relative momenta is strongly influenced by their 
mutual FSI and, for identical particles, also by QS. One can separate the FSI effect from 
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the production amplitude provided a sufficiently long two-particle interaction time in 
the final state as compared with the characteristic time of the production process. This 
condition requires the magnitude of the relative three-momentum q* = Pi — P2 = 2k* = 
Q in the two-particle cm. system much smaller than several hundreds MeV/c - the 
momentum transfer typical for particle production. For a two-particle bound state the 
momentum k* in this condition has to be substituted by (2/zeb) 1 / 2 , where e& is the binding 
energy and \i = m\m 2 j {m\ + m 2 ) is the reduced mass. 

Consider first the differential cross section for the production of a pair of non- 
identical particles 1 and 2 with the four-momenta Pi = {uji, pj} and the Lorentz factors 
7i = ujij nii. It can be expressed through the invariant production amplitudes T(pi, p 2 \ a) 
in the form 

(2^172 ,3^3 =El r (Pi.m°)l 2 . (2) 
d J pid J p 2 a 

where the sum is done over a = {S, M, a'}, i.e. the total spin S of the pair and its 
projection M (which is equivalent to the sum over helicities of the two particles) and 
the quantum numbers a' of other produced particles, including integration over their 
momenta with the energy-momentum conservation taken into account. 

We are interested in the pairs (1, 2) of the particles produced with a small relative 
velocity in a process with an ordinary phase space density of final-state particles so that 
a main contribution to the double inclusive cross section comes from the configurations 
(1,2, ..) with large relative velocities of the particles 1 and 2 with respect to other 
produced particles (i = 3, 4, ..n). Due to a sharp fall of FSI with the increasing relative 
velocity, we can then neglect the effect of FSI in all pairs and (2,i) except (1,2) 
and, in accordance with the upper diagram in figure [lj write the production amplitude 
as (see however [23] and subsection 4.3.3 for the account of the residual Coulomb field) 

T(p 1 ,p 2 ; a) = T (p 1 ,p 2 ; a) + AT(p 1 ,p 2 ; a). (3) 

Here To(pi,p 2 ;a) is the production amplitude in the case of no FSI, and AT(pi,p 2 ; a) 
represents the contribution of the FSI between particles 1 and 2, described by the 
formula 



AT{p 1 ,p 2 ;a) = — 1 - dK— — — j— — — , 4 

27r d J [k 2 — — iO)[(F — k) z — m 2 l — 1OJ 

where P = 2p = pi + p 2 , T (k,P — K]a) is the production amplitude analytically 
continued off mass-shell, f s {pi,P2] k,P — k) is the scattering amplitude of particles 1 
and 2 also analytically continued to the unphysical region. In the case of small k*, we 
are interested in, the central forces dominate so the scattering amplitude f s is diagonal 
with respect to the total spin S and doesn't depend on its projections. Since most of the 
systems of our interest (tt + tt~ , iC ± 7r T , 7i~p, K + K~ , K~p, except for pp) is described by 
a single value of S, we will often skip it to simplify the notation. 
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Let us express the amplitude T in a form of the Fourier integral 

T ( Pl ,p 2 ;a)= f d i x 1 d A x 2 e- ip ^- ip ^T(x u x 2 ;a) 

r J (5) 

= / d 4 xe-^ x/2 r P (x;a), 

where the last expression arises after the integration over the pair cm. four-coordinate 
X = [{p\P)xi + (p 2 P)x 2 }/ P 2 (d 4 xid 4 a;2 = d 4 Xd 4 a;) based on the separation of the 
phase factors due to the free cm. and relative motions: e~ ipiXl ~ mX2 = e - lPX e -^ x / 2 , 
Here the relative four-coordinate x = {t, r} = x\ — x 2 and the generalized relative four- 
momentum q = q — P(qP)/P 2 , q = P\ — p 2 ; note that qP = — m 2 2 . Apparently, 
the function T(x\, x 2 \ a) represents the production amplitude of particles 1 and 2 at the 
space-time points x\ and x 2 , respectively. It should be stressed that the representation 
(jSJ) concerns virtual particles as well. Inserting now in (j3J) the representation (j3J) with 
the substitutions p\ — > k, p 2 — ■> P — k, we get 

r( Pl , J , 2 ;a) = /d^d^*«-)(x 1 ,x 2 )r (ll ,x 2 ; Q ) 
d 4 xijj~^\x)Tp(x; a), 



where 



*$£(*i>*2) = *2£(^.^) " = e iPX ^| (+) (x) ' (7) 



coincides with the Bethe-Salpeter amplitude in continuous spectrum [21]. The second 
equality in (jSJ), similar to the one in (jSJ), merely arises after the integration over the 
pair cm. coordinate X as a consequence of the factorization of the free cm. motion in 
the phase factor e~ lPX . Thus, on the assumption of the quasi-free propagation of the 
low-mass two-particle system, the momentum dependence of the two-particle amplitude 
is determined by the convolution of the reduced production amplitude 

rp(x-a) = f d 4 Xe-' lPX T(x 1 ,x 2 ;a) (8) 

and the reduced Bethe-Salpeter amplitude ip~^~\x), the latter depending only on the 
relative four-coordinate x and the generalized relative four-momentum q. Using AS])-®, 
we can write 



where the correction to the plane wave is 



ipZ w (x) = e^ /2 + Atfj~ (+) (x), (9) 



i y ' 2vr 3 i 

,4 e iKX f s (p 1} p 2 ;K,P - k) 

K (k 2 - mi 2 + iO)[(P - k) 2 - m 2 2 + iO]' 1 ' 

In the two-particle cm. system, where P = 0, q = {0, 2k*}, x = {t*, r*}, the amplitude 
ip-^ (x) at t* = t* —t 2 = coincides with a stationary solution ^-k* ( r *) of the scattering 
problem having at large r* = |r*| the asymptotic form of a superposition of the plane 
and outgoing spherical waves [25J. 
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We see that one and the same production amplitude T(xi,x 2 ; a) or r P (x;a), 
corresponding to the space-time representation, enters into relations (jSJ) and ([6]). The 
effect of FSI manifests itself in the fact that the role of the functional basis, which the 
asymptotic two-particle state is projected on, is transferred from the plane waves to 
Bethe-Salpeter amplitudes ^>^~l(xi, x 2 ) or \x). 

Equation (jSJ) is valid also for identical particles 1 and 2 provided the substitution 
of the non-symmetrized Bethe-Salpeter amplitudes (xi, x 2 ) by their properly 

symmetrized combinations satisfying the requirements of QS: 

^-l(x 1 ,x 2 ) -> ± [tfJH^xa) + (-l)^-)^!,^)] . (11) 

In this case mi = m 2 , q — q and X = (xi + x 2 )/2. Similar to the case of non-identical 
particles, the assumption of the ordinary phase space density of the final-state particles 
allows one to account for the FSI and QS effects only in a given pair of identical particles 
produced with a small relative velocity. 

After substituting the representation (5) into (j2J), the double inclusive cross section 
takes on the form 

27r 6 7l72 

d 3 pid 3 p 2 

= J2 I d * x i d4x 2 dM dV 2 p P5 (xi, x 2 ; xi, x' 2 )^n x 2 )^H*K, x' 2 ) 
s J 

= W d 4 xd 4 x'p P5 (x;x')^ ( - ) (x)^ ( - ) *(x / ), (12) 
s J 

where the functions 

p PS (x 1 , x 2 ; x[, x' 2 ) = T{ x u x 2\ S, M, a')T*(xi, x' 2 \ S, M, a') 

M,a' 

p PS (x; x') = J2 M x \ S, M, a')r P (x'; S, M, a') = J d 4 X d 4 X' ( 13 ) 



M,a 



™W S (X + ^x, A" - A" + X' - Elf 



X 



p2 

represent elements of the unnormalized two-particle space-time density matrices; the 
density matrix p PS {xi, x 2 ; x[, x' 2 ) depends on the pair four-momentum P due to the 
account of the energy-momentum conservation in the sum J2 a '- O n the assumption 
of an instantaneous emission in the two-particle cm. system (t* = t 2 ), the second 
expression in f[T2"j) reduces to the ansatz used in [TUl [H] . 

Switching off the FSI and QS effects, for example, by mixing particles from different 
events with similar global characteristics, one can define the reference differential cross 
section 

fo x6 d6(T 

(2^)^x72 

d^pid rf p 2 

= E / ^ d M d 4 x 2 p P5 (xi, x 2 - xi, x^e-^^ 1 -^-^^-^) 



5 



X! / d A x d*x'p PS (x; x^^ )/2 (14) 
s J 
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and rewrite (TT21 as 

dV 
d 3 p!d 3 p 2 



E^(pi,p 2 )<^(^^ 2 )^K^2)>L, s (15) 



d 3 Pid 3 p 2 Y X P1 ' P P1 ' P iy /pip 2 5 



d 3 pid 3 p2 Efc(pi,P 2 )(^-(^^ ' ^ y/?ps 

where we have introduced the quasi-averages of the bilinear products of the Bethe- 
Salpeter amplitudes: 



P1P2S 



_ J d 4 a;i d 4 x 2 dV x d 4 4p PS (gi, x 2 ; ^ a^ggj (x 1; x 2 )^g-)*K, x 2 ) 

Jd 4 xi d 4 x 2 dV x d 4 4p P5 (xi, x 2 ; a^)e _ipi ^ 1 -<) -^2(^-4) (16) 

= U 8 ^ (xU s ^* (x>) )' = S^^PpsW)^-^)^^) 
\% 1 >/qP S J d*xd 4 x> 'pps(x;x')e-^ x - x 'y 2 

and the statistical factors Qs - the population probabilities of the pair spin-S* states in 
the absence of the correlation effect: 

Gs(Pi,P2) = 

J d 4 xi d 4 x 2 d 4 x[ dV 2 p PS (xi, x 2 ; x[, x ' 2 ) e ^p^-<)-^2-x' 2 ) 



EJ d%i d 4 x 2 d 4 xi d 4 x' 2 p PS (x u x 2 ; x[, x ' 2 ) e -w^-x\)-i P2 {x 2 -x' 2 ) 
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(17) 



J d 4 x d 4 x'p PS {x; x ') e -^-x')/2 
~ £ / d 4 x dVp PS (a;; xOe-^-*')^ ' 

Note that for unpolarized particles with spins ji and j 2 one has 

£ s = (2S + l)[(2i! + l)(2j 2 + 1)]- 1 = 1. (18) 

Generally, the spin factors are sensitive to particle polarization. For example, if two 
spin-1/2 particles were emitted independently with the polarizations Pi and P 2 , then 
Go = (1 - Px • P 2 )/4 and Q x = (3 + V x ■ P 2 )/4. 

The same procedure can be also applied to describe the production of weakly bound 
two-particle systems, like deuterons or hadronic atoms (tt + tt~ atoms, in particular). Due 
to a low binding energy, as compared with the energy transfers at the initial stage of 
the collision, there is practically no direct production of such bound systems. Their 
dominant production mechanism is thus due to the particle interaction in the final 
state. The invariant production amplitude Tb(Pb, S, M, a') of a spin-S" bound system 
b = {1 + 2} is then described by the lower diagram in figure [1] corresponding to the 
second term in the upper diagram with the free two-particle final state substituted 
by the bound one. Therefore, similar to (JSJ), this amplitude is related to the Fourier 
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transforms T or Tp b of the off-mass-shell two-particle amplitude T : 

T b (P b ; S, M, a') = f d 4 x x d 4 x 2 ^ { pJ(x 1 ,x 2 )T(x 1 , x 2 ; S, M, a') 

r ' (19) 

= j d 4 x^ { -\x)T Pb (x;S,M,a') 

and the corresponding differential cross section - to the same two-particle space-time 
density matrices p PS or pps as enter into (EES]) , up to the substitution P —> P b : 



(2tt)S 



3„, u u b 



d 3 P b 

| d 4 X! d 4 X 2 dM d%p PbS (x 1 , X 2 - X[, x' 2 )^ Pb \ Xl , X 2 )<y Vl, 4) ( 20 ) 

d 4 x dVp n5 (a;; x')ip b { ~' (x)i/) b { ~>*(x'). 



Here ^^(xi, x 2 ) = [^fp/ (^l, ^2)]* = [e^^^f (2)]* is the Bethe-Salpeter amplitude 
for the bound system. At equal emission times of the two particles in their cm. system, 
the amplitude ip b describing their relative motion, coincides with the usual non- 

relativistic wave function in discrete spectrum ?/>f (r*). Similar to ( fl5i) . one can also 
rewrite the production cross section of the bound system through the reference cross 
section in (TH1) taken at pi = V b m\j (mi + m 2 ), p 2 = V b m 2 j(rn\ + m 2 ) and 71 = 7 2 = 76 
(i.e. q = 0, P = P b ): 

dV UA) 



'd 3 Pl d 3 P2 v ^'^ \ Tb v /r " ^'/ops' 
We see that the production of a weakly bound system {1 + 2} is closely related with 
the production of particles 1 and 2 in continuous spectrum at small kinetic energies 
in their cm. system. This relation was first formulated [26] in connection with the 
production of non-relativistic deuterons and then generalized [27J to the relativistic 
case and the inclusive production. Similar relation was obtained, in the limit of an 
instantaneous emission from a point-like region, also for the case of the production of 
pure Coulomb hadronic atoms [10] . A complete treatment of the production of weakly 
bound systems, accounting for the finite-size effect, can be found in [TT] . 

3. Approximate description of the 7r + 7r production 



Following [10] , let us first neglect the ir + n~ strong FSI and assume only two types 
of pion emitters: SL emitters (e.g., p- or A-resonances) characterized by small sizes 
or decay lengths on a fm level, and LL emitters (e.g., rj, K s or A) with very 
large or macroscopic decay lengths. Since the relative space-time distance between 
the emission points x enters in the pure Coulomb amplitudes ^~) coul (x) and 

^(-)coui^ gca ] ec j by ^ e B orir radius a = —387.5 fm, one can put in ([12]) and (I20I) 
^ _) (x) « e-^^-^^O) = e ik * r *[^i(0)]* and ip { b ~\x) « ^~ )c ° ul (0) = [^ coul ( )]* 
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for the fraction A of the pairs with both pions from SL emitters (r* <C \a\), and 
^(-)coui^ ^ exp[ik*r*— i(/c*a)~ 1 ln(k*r* + A;*r*)] (a plane wave amplitude with the phase 
modified by Coulomb interaction) and ipl ^ coul (a;) 



^- )C ° Ul (oo) 



for the remaining 



fraction (1 — A) of the pairs with at least one pion from a LL emitter (r* ^> \a\). As a 
result, and (J2D5 reduce to: 

d 6 a d 6 d 



d 3 pid 3 p 2 

dV 6 



d 3 P fc 



Wit 



d 3 pid 3 p 2 
dV 



A 



i(0) +(1-A) 



d 3 pid 3 p 2 



A 



^ COUl 



(0) 



(22) 



(23) 



where <To represents the production cross section of the non-interacting pions and the 
expression for the production of bound 7r + 7r~ system implies Pi = p 2 = Pft/2 and 
7i = 72 = 76- The squares of the non-relativistic Coulomb wave functions at zero 
separation are well known: 



7] = (k*a) 



-i 



A coul (o) 



5. 
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(7r|a| 3 n 3 ) 



{nl}, 



(24) 
(25) 



where the Bohr radius a = (/xe^) -1 is negative for 7r + 7r~ system due to the opposite 
signs of 7r + and tt~ charges (ei = — e 2 = e). The Coulomb penetration factor A c {j]) 
(sometimes called Gamow factor) behaves at small Q = 2k* < Q c as Q c /Q and at large 
Q approaches unity as 1 + \Q C /Q, where Q c = 4ir/\a\ = 6.4 MeV/c. As for the bound 
7r + 7r~ states b = {nl}, only the s-wave states {n0} are produced at zero separation and 
their fractions with given main quantum numbers n are uniquely fixed by the n~ 3 law 
in (EH]). 

The numbers Na and iVjf of produced and breakup 7t + tt~ atoms, required to 
calculate the breakup probability ([T|), can then be obtained in two steps [ID]. First, 
one simulates the non-correlated two-pion spectrum d 6 -/Vo/d 3 pid 3 p 2 or, constructs it by 
mixing pions from different events, and determines the overall normalization parameter 
g and the fraction A or A = Xg by fitting the theoretical spectrum 

d 6 N d 6 iVo 



d 3 pid 3 p 2 

d 6 iVo 
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d 3 pid 3 p 2 
[AA c {n) + A' 



[XA C ( V ) + (1 - A)] 



(26) 



d 3 pid 3 p 2 

to the measured spectrum of the pion pairs; to get rid of the pairs from the breakup 
of the 7r + 7T _ atoms in the target, the fit should be done in the region Q >~ 3 MeV/c. 
In the second step, one can use ( l23l) and the fitted parameter A = Xg to calculate 
the three-momentum distribution of the numbers of produced atoms in given states 
b = {n0}: 



d 3 N b 
d 3 P„ 



(27r) 3 7b 



d 3 iVn 



d 3 pid 3 p 2 



A 



n°0 Ul (0) 



(27) 
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where Pi = P2 = P&/2. Then, taking into account that d 3 pid 3 p 2 — ^d 3 Pd 3 q = 
^•yd 3 Pd 3 Q and that at sufficiently small Q < Qo the non-correlated two-pion spectrum 
is practically Q-independent (up to a correction 0{Ql/rrv^)), one can calculate Na from 
the number Nq(Q < Qo) of simulated or mixed non-correlated pion pairs with Q < Qq\ 

One can also calculate Na from the number of correlated non-atomic pairs, 

KT n -(Q < Qo) = ^N (Q < Q ) f Q ° dQQ 2 A c ( v ) } (29) 
using so called fc-factor [4]: 

HQ < Qo ) = jw^S- (Q < <M - ig,^ (30) 

e.g., fc = 0.615, 0.263 and 0.140 for Q = 2, 3 and 4 MeV/c respectively. As for the 
number of breakup atoms A^jf , it is simply obtained by subtracting the fitted numbers 
of correlated (cna) and non-correlated (nc) non-atomic pion pairs (see (1261) ) from the 
measured number of pion pairs: 

N% = AW (Q < Q cut ) - AX- (Q < Qcut) - (Q < Qcut) (31) 

A£V(Q < Q cut ) = A'N (Q < Qcut). (32) 

The value of Q cnt should be chosen sufficiently large so that the interval (0, Qcut) contains 
the signal from practically all atomic pairs. The possible choice is Q cut = Qo = 4 MeV/ c. 
One can also choose a smaller value, correcting for the loss of atomic pairs with the help 
of simulated efficiency factor e^ r [I] . One can obtain A^jf also in a more direct way using 
the data from multi-layer targets [2J. 

Let us now consider the modification of (!26j) and (j27|) due to the strong FSI and 
finite space-time separation of the particle emitters. Formally one can write 



d 6 N . d 6 A r 



{A[l + 5(k*)]A c (77) + A'} (33) 



d 3 pid 3 p 2 d 3 pid 3 p 2 

(34) 

iV.^iV (Q<Q )(^-) 3 E^(l + ^) (35) 
KZ-(Q < Qo) = ^N (Q < Qo) [ Q ° dQQ 2 A c ( V )[l + S(k*)} (36) 

Qn JO 



Ql 

(47T) 3 En"- 3 (1 + ^) 

47r 2 |a| 3 / Qo dQ Q 2 A c (r))[l + S{k*)} ' 



*« < <w - w:.: r:" r , (37) 
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1 + 6Qf) = (vf^^V))-^ IMvT 1 (38) 
1 + 5 n = (^\x)^>(x'))'* p |C°o ul (0)f 2 • (39) 



where the correction factors are determined by the averaging of the bilinear products of 
the reduced Bethe-Salpeter amplitudes over the distribution of the relative space-time 
separations of the SL emitters: 

'SL 
qP 

r':V,' u (0) 

OP 

The averaging is defined in ( TT6l) with the reduced space-time density matrix substituted 
by its part, p S p , related only with the SL emitters. Equations (138]) and (1391) account 
only for the elastic transition a — > a and ignore a small contribution of the inelastic one 
f3 — > a, where a = {/T + 7r~}, (3 = {7r°7r }; see section 6 and (11611) . (I162p for the complete 
treatment. 

In fact, it can be argued [13J, [H] that 

5 n « 5(0) (40) 

provided the characteristic spatial separation of the pion SL emitters in the two-pion cm. 
system is much less than the two-pion Bohr radius \a\. This result immediately follows 
from the well known Migdal's argument [26]. Namely, since the particles in continuous 
spectrum at zero kinetic energy and in discrete spectrum at very small binding energy 
k 2 /{2p) — > are described by practically the same wave equations, the r*-dependence 
of the corresponding wave functions at a given orbital angular momentum should be the 
same for the distances r* <C k _1 (i.e. r* <C n\a\ in the case of a hadronic atom with the 
main quantum number n). 

One may see that the approximate equality ( HOI , together with the assumption 
of a weak /c*-dependence of the correction <5(k*), justify the use of the approximate 
equations ( |26l) - (l30l) . In this approximation, the finite-size correction merely reduces to 
the rescaling A -> A[l + 5(0)]. 

We show in sections 5 and 6 that ( HOI is subject to a normalization correction 
0(/o/ a ) ~ 0.3%/n and other small corrections 47rO(fc| ( (/ 3Q! ) 2 /a) and 0(a~ 2 ), where 
kp = 35.5 MeV/c is the momentum in the channel j3 = {7r°7r } at the threshold transition 
to the channel a = {7r + 7r~}, f$ a = y/2(ciQ — a[])/3 ~ —0.2 fm is the transition amplitude. 
Taking the normalization correction explicitly into account (see (11541) . (11431) and (I147P ): 

1 + 5 n = (1 + 5' n ) 



\ +<>(n) Vn 



0((A£/, 



/3c^ 2 



n M (41) 

where <fi(n) m 3 is defined in ( 1109[) . one can rewrite the approximate equality in (|4~0~|) as 

5' n = 5(0) + 47rO(k;(C) 2 /a) + 0«r* 2 /a 2 > SL ). (42) 

The neglect of the corrections 5 n and 5(k*), i.e. the use of the approximate equations 
(|26|) - (|30|) instead of (!33|) - (l37|) . leads to the systematic shift of the breakup probability: 

AP br /P br = -AN A /N A + AN h //N h A \ (43) 
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To estimate this shift, one can approximate the correlation function 

where S(k*) = (5(k*)) is the finite-size correction averaged over the pion three-momenta 
at a fixed k* = |k*| = Q/2, by 

K{Q) = AA C ( V ) + A' (45) 

and use the fitted parameters A, A' to calculate ANa/Na and AA^/A^f: 

AiV A Enn -3(i + 5n _A) 



N A E n n~ 3 (l + S n ) 

an% _ n^AQ < gcuO J Qcut dgg 2 [^(g) - g(g)] 
a^ A^(g<g cut ) ; Qcut dgg 2 ^(g) 



(46) 



(47) 



In the experiment DIRAC, N%? v -/N% = 16 at Q cnt = 4 MeV/c [4j and, for Q cut < Q c , 
when QA c {rj) « const, it decreases approximately quadratically with decreasing g cut 
provided the signal region contains most of the atomic pion pairs, i.e. down to Q C ut ~ 2 
MeV/c. 

In the following, we perform an analytical and numerical study of the corrections 
5 n and S(k*) and their effect on the breakup probability Pbr- Here we only mention that 
the condition r* <C n\a\ can be violated for pion pairs containing pions from the decays 
of some resonances such as u and rj' with the decay lengths of about 30 and 900 fm 
respectively. The corresponding exponential tails are clearly seen in figure [2j where the 
redistribution simulated with the UrQMD transport code [28] is shown for pion pairs 
produced in pNi interactions at 24 GeV in the conditions of the DIRAC experiment at 
CERN [291. 



4. Correlation femtoscopy — basic assumptions 

4-1. Non-interacting non-identical particles: space-time coherence 

To clarify the meaning of the two-particle space-time density matrix p s (x\, x 2 ; x[, x' 2 ), 
let us first neglect the FSI effect and substitute the Bethe-Salpeter amplitudes by the 
plane waves. Changing in (fl2"|) the integration variables Xi,x\ by the new ones: 

Xi = -(xi + x'i) 6i = Xi- x' p (48) 

we can rewrite the production cross section of two non-identical particles as 

d% 



(27r) b 7l 72 



d 3 P id 3 p 2 (4g) 



J d 4 xi d A x 2 G s (xi, pi, x 2 ,p 2 ) = 12 J d 4 xg PS (x, q), 
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where X — X\ — X2 and the real functions 

G s (xi,p 1 ;x 1 ,p 2 ) = 

J d 4 ei d 4 e 2 e- ipiei - ip2e2 p P5 (xi + |,^+ - - 

( ~\ r^irn {y^P* P - t P lP - \ (50) 
9ps{x; q) = J d IG S IX + -j^x, Pi;X - -j^x, p 2 j 

The function G5, usually called emission function, being a partial Fourier transform 
of the space-time density matrix, is closely related to the Wigner density, the latter 
collecting all contributions due to free streaming of the emitted particles to given space- 
time points through an integral over the emission function (see (49) in |30j). 

It is clear from fj49l) and fl50l) that more narrow is the width of the diagonal of the 
space-time density matrix (the width of the ej-distribution), more wide is the distribution 
of particle four-momenta. In particular, the diagonal space-time density matrix (i.e., 
zero width of the ej-distribution) would yield the uniform four-momentum distribution, 
in correspondence with the infinite uncertainty in the four-momenta of the particles 
localized at certain space-time points. 

Consider as an example the particle emission by independent one-particle emitters 
of various types A according to the one-particle production amplitudes (see also [3T] ) 

(xi - x A ) 2 (x i - x 0A ) 2 



T^\xi, x A ) ~ v A (x A ) exp 



v A (x A ) ~ exp. ., ., 
\ 4r n 4rr / 




(51) 



These amplitudes correspond to the emitters at rest with a Gaussian distribution of the 
emission points X\ = {t±, i^} around the emitter centers x A = {t A , r A }, also distributed 
according to a Gaussian law. In four-momentum representation, 

Tx\pi] x A ) ~ v A (x A )u A (pi) exp(-zpix A ) 

u A (pi) ~ exp(-r A 2 pi 2 / 2 - ta 2 Poi 2 / 2 ) • 

Assuming further that the emitters are sufficiently heavy, we can describe them 
classically. The four-coordinates of the emitter centers x A can then be considered as a 
part of the quantum numbers a', the sum in (1131) thus containing the integration over 
x A . Performing this integration, we get for the elements of the one-particle space-time 
density matrix related to the emitter A: 



p ( A l \x 1 ,x' 1 ) = J d 4 XA / 7 A 1) (xi;x y i)'7 A 1) *(x / 1 ;a;A) 



ex 2 e 01 2 \ / x x 2 x 01 2 \ ( 53 ) 



ex P -7—5- 7-^ ex P 



4r A 2 4r A 2 / *\ 2r 2 + r A 2 2t 2 + t a 2 
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and for the corresponding emission function: 



G { }\x 1 ,p 1 ) = J d 4 e ie -^V« (xt + 



£1 
2 



) 



(54) 




The contribution of the emitter A to the single-particle production cross section is 



We may see that the emitter space-time dimensions and ta determine both the space- 
time coherence of particle production (the non-diagonality of the space-time density 
matrix) and the distribution of particle four-momenta. In particular case of the emitters 
of a vanishing space-time extent: r A = TA = (no coherence) any particle four- momenta 
are equally probable. 

Note that for the emitter moving with a non-relativistic velocity (3 A and emitting 
a particle 1 with the mean three- momentum = mif3 A , the amplitude fl5T|) and 
the density matrix ( 1531) respectively acquire phase factors e _ipAXl and e~ ipA€l and the 
substitution pi — > pi — has to be done in the amplitude ua(pi)- After averaging over 
the p^-distribution that decouples from the distribution of other emitter characteristics 
in a Gaussian form of a width A , we still arrive at fl53|) - fl55|) . up to a substitution 
r a' — > TA 2 /[2(rAAo) 2 + 1] in the e- and momentum-dependent factors, corresponding to 
a widening of the momentum distribution due to the dispersion of the emitter velocities. 

As for the actual values of the parameters and ta, we can estimate them using the 
information about particle transverse momenta, pt, which are much less influenced by 
the motion of the emitters than the longitudinal ones. Doing this for pions or kaons, we 
should however exclude the \ow-p t region which is dominantly populated by the decays 
of low-lying resonances. We can also use the ^-distributions of these resonances. In 
both cases the pt 2 -slopes in the interactions of elementary hadrons are of ~3 (GeV/c)~ 2 
(see, e.g., [32]), yielding on average ta 2 ,T4 2 ~ 0.1 fm 2 . Somewhat larger values can 
be expected in heavy ion collisions where a substantial part of the emitters can be 
associated with the centers of the last rescatterings characterized by sufficiently large 
momentum transfer. It is important that the estimated values of r^ 2 , T4 2 appear to be 
much smaller than the effective values of the parameters ro 2 , tq 2 obtained in femtoscopic 
measurements. The latter being of about 1 fm 2 for pions produced at p t ~ (pt) in 
hadron-hadron interactions and up to several tens fm 2 in the collisions involving heavy 
nuclei. 

4-2. Non-interacting identical particles: QS correlations 

4-2.1. Correlation function Consider the production of non- interacting identical 
particles. It should be noted that this consideration is not of academic interest only. 
Thus for identical pions or kaons, the effect of the strong FSI is usually small and 




(55) 
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the effect of the Coulomb FSI can be in first approximation simply corrected for (see 
[33] and references therein). The corrected correlation effect is then determined by the 
QS symmetrization only i.e. the Bethe-Salpeter amplitudes have to be substituted by 
properly symmetrized combinations of the plane waves (see (fTTT) ). As a result of the 
interference of these waves, there appears the additional term, not present in (1491) : 

d 6 rr 
d' i pid J p 2 

= E /d^ 1 d^[ Gs (, 1 ,p 1 ;, 2 ,p 2 ) + G s (x 1 , K x 2 ,p)(- 1 )«oos W ] (56 ) 
s J 

Note that the off- mass-shell four- momentum p = (pi + p 2 )/2 enters as an argument of 
the emission function G$ in the interference term. 

It is convenient to define the correlation function lZ(pi,p 2 ) as the ratio of the double 
inclusive cross section d 6 <r to the reference one d 6 o"o which would be observed in the case 
of absent QS and FSI effects: 

K(p " P2) = d%te^)' (57) 

In the high-energy collisions involving nuclei, we can neglect the kinematic constraints 
as well as rather weak dynamical correlations and construct the reference distribution 
using the particles from different events with similar topology. In case of a negligible 
FSI, there is no correlation for non-identical particles: 7t{pi,P2) — 1> while for identical 
particles the correlation arises due to the interference effect: 

J2 J d 4 xi d 4 x 2 G s (x l} p; x 2 ,p)(-l) s cos(gx) 

Kpi,Pz 1+ T,fd 4 x 1 d 4 x 2 Gs(x 1 ,p 1 ;x 2 ,P2) 

s 

= i + Y.Gs(-i) s (cos(qx)y; iP2S 

S - ( 58 ) 

E / d 4 xg PS (x, 0)(-l) 5 cos(gx) 

= 1 + ^ 

£/ d 4 xg PS (x,q) 
s 

= l+E^(-l) 5 (cos(ga;));' P5 , 
s 

where the quasi-averages satisfy the equalities 

{cos{qx)) pip2S -{e } ^ 

= (cos(gx));' P5 =(e^')/ 2 )^; 

the factors Q$ represent the population probabilities of the pair spin-S* states out of the 
region of the correlation effect. They are defined in (|17p and can be expressed through 
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the emission functions as 

/ d 4 xt d 4 x 2 G s {x 1 ,pi; x 2 ,p 2 ) 



Gs(Pi,P2 



J2s I d 4 xi d 4 x 2 G s (x 1 ,p 1 ; x 2 , p 2 ) 
Jd 4 xg PS (x,q) 



4_ r„ (60) 



EsJd 4 xg PS {x,q) 

They can be also considered as the initial (QS switched off) statistical factors. For 
initially unpolarized spin-j particles: J2sGs(— l) 5 = ( — l) 2j / (2j + 1). 

Assuming, for example, that for a (generally momentum dependent) fraction A of 
the pairs the particles are emitted by independent SL one-particle emitters described by 
the Gaussian amplitudes (I5TT) or (1521) . while for the remaining fraction (1 — A), related 
to LL emitters (rj, K®, A, . . . ), the relative distances r* between the emission points in 
the pair cm. system are extremely large, the correlation function 

^(Pl>P2) 

, L , W/ 1 ,5 ^EA,B^(pi)^(p 2 )^(p 2 )^(p l )e-^-^) (61) 
s Ea,b\ua{Pi)u b {P2)\ 2 
where the sum J2a,b is done over all characteristics of the emitters. In the case of only 
one type of the SL emitters that are at rest and differ only by the four-coordinates xa of 
their centers, the amplitudes ua reduce to a single universal amplitude u and the sum 
merely reduces to the averaging over xa, i.e. 

K(p u p 2 ) = l + \Y,Gs(-l) S (cos(q(xA-x B ))) 

S 

= 1 + A£fc(-l) 5 exp(-r V - r 2 g 2 ) . ^ 
s 

We see that a characteristic feature of the correlation function of identical particles is the 
presence of an interference maximum or minimum at small |q|, changing to a horizontal 
plateau at sufficiently large |q|, large compared with the inverse characteristic space-time 
distance between the particle emission points. 



4-2.2. Smoothness assumption In the simple model of only one type of the emitters 
contributing to the observable interference effect and in the absence of the relative 
emitter motion, the width of the low-|q| structure is solely determined by the 
characteristic space-time distance between the one-particle emitters and does not depend 
on the parameters and ta, characterizing the space-time extent of the emitters 
themselves - see ( 1621) . It means that the enlargement of the production region 
related to the latter (r 2 — > r 2 + |^a 2 , t o 2 —> T o 2 + \ t a 2 ) is compensated in the 
correlation function due to the different momentum arguments of the emission functions 
in the numerator and denominator of fl58|) . This is clearly seen when calculating the 
correlation function directly from (J2J), substituting the production amplitude T(pi,p 2 ; a) 
by the symmetrized product of the Kopylov-Podgoretsky one-particle amplitudes in 
momentum representation, see (|6TT) and ( 1621) . Of course, the independence of the 
interference effect on the space-time extent of the emitters in this model (assuming 
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that the emitters decay according to a single universal amplitude u and differ by the 
four-coordinates of their centers only) is justified only in the case of sufficiently small 
overlap of the emitters to guarantee the assumption of their independence. 

Generally, even in the case of independent emitters, the particles are emitted 
by moving emitters of different types and the correlation function depends also on 
their space-time extent r^, ta- Particularly, for a Gaussian distribution of the mean 
emission three-momentum of a width A , f l62|) is modified by the substitution [31] 
r 2 — > r 2 + rA 2 /[2 + (r^Ao) -2 ]. Usually, the effect of a finite space-time extent of the 
one-particle emitters is negligible: 



Note that these conditions guarantee sufficiently smooth four-momentum dependence of 
the emission function Gs(xi,Pi',X2,P2), such that we can neglect its dependence on the 
four-momentum difference q in the region of the interference effect characterized by the 
inverse space-time distance between the particle production points. On this, so called 
smoothness assumption, ( J58|) reduces to: 



Equation (JM]) is valid up to a correction representing a fraction of ryi 2 /r 2 , ta 2 /tq 2 . 
This correction composes a few percent for high-energy hadron-hadron collisions and a 
fraction of percent for the collisions involving heavy nuclei. Note that fl64|) is often used 
to calculate the correlation functions of non-interacting identical particles with the help 
of various classical transport codes (like RQMD, VENUS or UrQMD) [2H] - the emission 
points are identified with the points of the last collisions or the resonance decays. 

At sufficiently small Q, one can calculate the one-dimensional correlation function 
7Z(Q) using a more simple and faster procedure than the averaging according to (164|) . 
For this one can exploit the fact that the angular distribution of the vector Q becomes 
isotropic at Q — > and calculate (cos(gx)) = (cos(Qr*)) by averaging over the uniform 
distribution of the cosine of the angle between the vectors Q and r* and, over the one- 
dimensional redistribution determined at Q — > 0. We have checked the accuracy of 
this procedure using the UrQMD simulation of the pNi interactions at 24 GeV in the 
conditions of experiment DIRAC [2H] and determined the redistributions in various 
Q-intervals. It appears that the tt~tt~ correlation functions calculated from the r*- 
distributions corresponding to the intervals 50-100, 100-150 and 150-200 MeV/c agree 
with that corresponding to the redistribution in the lowest Q-interval of 0-50 MeV/c 
within 0.2-0.9%, ~ 3% and ~ 7%, respectively. It should be noted that the increasing 
difference of the correlation functions with the increasing lower boundary of the above 



r A 2 /2<r 2 , r//2<r 2 . 



(63) 



n( Pu p 2 ) = 1 + £(-l) 5 fe (cos(gx)) 



(64) 




(65) 
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Q-intervals is not related with the violation of the smoothness assumption but rather 
with the approximate treatment of the angular dependence of the vector Q and with 
the Q-dependence of the fractions of pairs containing pions from resonance decays. 

4-2.3. Femtoscopy with identical particles One can see from ( 1621) that, due to the 
on-shell constraint q Q = vq = vq^, the correlation function at vt q > r substantially 
depends on the direction of the vector q even in the case of spherically symmetric spatial 
form of the production region. Thus the transverse (q _L v) and longitudinal (q || v) 
correlation radii are Tt = ro and = {r 2 . + v 2 ^) 1 ^ 2 respectively. 

The on-shell constraint makes the g-dependence of the correlation function 
essentially three-dimensional (particularly, in pair cm. system, qx = — 2k*r*) and thus 
makes impossible the unique Fourier reconstruction of the space-time characteristics 
of the emission process. However, within realistic models, the directional and velocity 
dependence of the correlation function can be used to determine both the duration of 
the emission and the form of the emission region [15], as well as - to reveal the details 
of the production dynamics (such as collective flows; see, e.g., [341 [35] and reviews 
[TBI [191 1201 El]). For this, the correlation functions can be analyzed in terms of the 
out (x), side (y) and longitudinal (z) components of the relative momentum vector 
q — {<lx, q y , q z } [36], [37]; the out and side denote the transverse components of the 
vector q, the out direction is parallel to the transverse component of the pair three- 
momentum. The corresponding correlation widths are usually parameterized in terms 
of the Gaussian correlation (interferometry) radii r i; e.g., for spin-0 bosons 

n(p 1 ,p 2 ) = 1 + Aexp(-r^ - r 2 y q 2 - r\q 2 z - 2r 2 xz q x q z ), (66) 

and the radii dependence on pair rapidity and transverse momentum is studied. The 
correlation strength parameter A can differ from unity due to the contribution of 
LL emitters, particle misidentification and coherence effects. Equation ( l66l) assumes 
azimuthal symmetry of the production process. Generally, e.g., in case of the correlation 
analysis with respect to the reaction plane, all three cross terms q^qj contribute. 

It is well known that particle correlations at high energies usually measure only a 
small part of the space-time emission volume, being only slightly sensitive to its increase 
related to the fast longitudinal motion of particle emitters. In fact, due to limited 
emitter decay momenta p^ec of few hundred MeV/c, the correlated particles with nearby 
velocities are emitted by almost comoving emitters and so - at nearby space-time points. 
In other words, the maximal contribution of the relative motion to the correlation radii 
in the two-particle cm. system is limited by the moderate emitter decay length rp^ ec jm. 
The dynamical examples are resonances, colour strings or hydrodynamic expansion. To 
substantially eliminate the effect of the longitudinal motion, the correlations can be 
analyzed in terms of the invariant variable Q = 2k* = (— q 2 ) 1 ^ 2 and the components 
of the three-momentum difference in the pair cm. system (q* = Q = 2k*) or in the 
longitudinally comoving system (LCMS) [38]. In LCMS, each pair is emitted transverse 
to the reaction axis so that the generalized relative three-momentum q coincides with 
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q*, except for the cm£-component q x = 7^*, where 7 t is the LCMS Lorentz factor of 
the pair. Particularly, in the case of one-dimensional boost-invariant expansion, the 
longitudinal correlation radius in LCMS reads [35] 

r z = (T/m^r, (67) 

where T is the freezeout temperature, r is the proper freezeout time and m t is the 
transverse particle mass. In this model, the side radius measures the transverse radius 
of the system while the square of the out radius gets an additional contribution 
(pt/mt) 2 Ar 2 due to the finite emission duration At. The additional transverse 
expansion leads to a slight modification of the ^-dependence of the longitudinal radius 
and - to a noticeable decrease of the side radius and the spatial part of the out radius 
with p t . Thus in the case of a linear non-relativistic transverse flow velocity profile 
Pf — Po r t/R of the expanding fireball with the freeze-out transverse radius R, the side 
radius 

r y ^R/(l + mtPo/T) 1/2 . (68) 

The decrease of the two-pion correlation radii with increasing transverse mass 
(expansion) and decreasing centrality (geometry) has been demonstrated, e.g., in 
Au+Au collisions at i/s NN = 200 GeV [39] . 

Since the freeze-out temperature and the transverse flow determine also the shapes 
of the m t spectra, the simultaneous analysis of correlations and single particle spectra 
for various particle species allows one to disentangle all the freeze-out characteristics 
(see, e.g., [IE])- Thus in heavy ion collisions the correlation data show rather weak 
energy dependence and point to the kinetic freeze-out temperature somewhat below the 
pion mass, a strong transverse flow (with the mean transverse flow velocity of about 
half the velocity of light), a short evolution time of 8-10 fm/c and a very short emission 
duration of about 2-3 fm/c (see, e.g., a recent review [2~T]). 



4-3. Interacting particles: FSI correlations 

4-3.1. Production of interacting particles It is clear that the smoothness assumption 
allows one to express the production cross section through the emission function 
Gs{xi,Pi, x 2 ,p 2 ) also in the case of interacting particles. Thus, separating the two- 
particle cm. motion in the phase factor exp[iP(X — X')] = exp[z(pi — q/2)e\ + i(j) 2 + 
q/2)e 2 ] and using the smoothness assumption to neglect here q compared with p 12 lil and 

S(-\-) 

substitute, in the amplitudes ip~ (x), the relative coordinates x = x + (ei — e 2 )/2 and 
x' = x — (ei — 62)/ 2 by their mean value x, we can rewrite ( fl2|) in a simple approximate 
form: 



(2^) 6 7i7 2 /^- = E / d i x l d'x 2 Gs(x 1 , Pl ;x 2 ,p 2 ) ^ {+ \x) 
d^pid^ 3 J 



\ The account of q in the phase factor would lead to the substitution of the particle four-momenta in 
the emission function by their mean (off-mass-shell) values: pi — > Pmij{ra\ + TO2). 
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(2vr) 6 7l 72 



d 4 xg P s{x,q) ip~ i+ \x 
d% 



d^p 2 ^ gs 



^ + \x) 



qPS 



(69) 



where d 6 <7o is the production cross section of non-interacting particles introduced in 
019]). The averaging (. . -)qpg and the initial spin factors Qs are defined in ( 1651) and (I60p . 
The correlation function defined as the ratio d 6 o~/d 6 cr then takes on the form: 

2^ 



s 



(70) 



Recall that for identical particles the Bethe-Salpeter amplitudes ip~ (x) should be 
symmetrized according to (fTTT) . 

Note that for non-identical particles, one also arrives at (1691) and (1701) using the 
approximate ansatz ^(i^ij) = e 1 ^ 161 ^ 262 ^^^^!,^) which becomes exact in 
the absence of FSI. For identical particles, this ansatz, applied to the non-symmetrized 
amplitudes '$ rS ( + ), leads to the correlation function (see also (58) and (60) in [30] ) 
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where if) is the reduced non-symmetrized Bethe-Salpeter amplitude (V'q ( 
for non-interacting particles). Clearly, the smoothness assumption allows one to put 
(• • -YqPS — (• ■ -)qPS an d thus recover symmetrized equation (1701 . 

Similar to the case of non-interacting particles, the relative correction to the 
smoothness approximations in (l69i) -( 17Ti) is determined by the ratios r^ 2 / r Q 2 , ta 2 / t 2 - the 
measures of the non-diagonality of the space-time density matrix. For identical particles, 
the correction arises mainly from the simplified treatment of the symmetrization effect 
and, according to subsection 4.2.2, it is expected on a few per mil level for the processes 
involving heavy nuclei. For non-identical particles, the corrections to the finite-size FSI 
contributions are of the same order while, those to the complete correlation functions 
are usually substantially smaller, being scaled by the relative finite-size contributions 
of the strong and Coulomb FSI. In case of <C r* <C \a\, we are interested in, 
the corresponding strong and Coulomb FSI contributions are of 2f s /r* and 2r*/a 
respectively (see section 5). 

Proceeding in a similar way with the production cross section of a bound two- 
particle system, we arrive, on the same conditions as in the case of continuous spectrum, 
at the approximate form: 



(27t) 3 7 6 
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d 3 P* 



d 4 xi d A x 2 G s {xi,pi; x 2 ,p 2 

d 4 x5f P5 (x,0) ipb i+ \x 
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where p, = P b m i /(m 1 + m 2 ), P = P b and P = P b0 . 

4-3.2. Equal-time approximation For non-interacting particles, the non-symmetrized 
Bethe-Salpeter amplitude ip~~\x) = e _lk * r * is independent of the relative emission time 
t* in the pair cm. system. On the contrary, the amplitude of two interacting particles 
contains an explicit dependence on t* - the interaction effect vanishes at \t*\ — > oo. 
However, it can be shown [T2] (see Appendix A) that the effect of non-equal times can 
be neglected on condition 

\t*\ < m(t*)r* 2 , (73) 

where m(t* > 0) = m<i and m(t* < 0) = mi. On this condition one can use the 
approximation of equal emission times of the two particles in their cm. system (t* = 0) 
and substitute the Bethe-Salpeter amplitude by the usual non-relativistic two-particle 
wave function. The applicability condition (1731 of the equal-time approximation is 
usually satisfied for heavy particles like kaons or nucleons. But even for pions this 
approximation merely leads to a slight overestimation (typically less than a few percent) 
of the strong FSI contribution to the production cross section [12]. To demonstrate 
this, one can use the simple static Gaussian model of independent one-particle emitters 
described by the amplitude (j5ip . The applicability condition (1731 of the equal-time 
approximation can then be written as [T2] 

r < M7?"o(ro 2 + -A) 2 ) 1/2 - (74) 

Recall, however, that in high-energy collisions, the static model is relevant for a limited 
rapidity region only. It means that the pair velocity v in the rest frame of the 
contributing emitters is essentially determined by the distribution of particle transverse 
momenta. For pion pairs at Q — > one then has (v) ~ 0.8. For r <~ r condition 
( 1741) requires sufficiently small Compton wave lengths of the particles in the emitter rest 
frame: 1/cUj <C r , while for large characteristic emission times, r ^> r Q /v, it requires 
small de Broglie wave lengths: 1/pi <C tq. Clearly, this condition is not satisfied for very 
slow particles emitted by the emitters of a long lifetime. The increasing importance of 
the non-equal time effect with the decreasing pair velocity and increasing lifetime of 
the emitters is demonstrated in figures [3] and H] for the FSI contribution in the 7r°7r° 
correlation function. For sufficiently large velocities v > 0.5 and radii r Q > 1 fm, we are 
interested in, the effect is rather small, not exceeding 5% of the FSI contribution in the 
low-fc* region, corresponding to the effect of a few per mil in the correlation function. 

As for the effect of non-equal times on the Coulomb FSI it doesn't influence the 
leading zero-distance (r* = 0) part and, the effect of the subleading part (expected on a 
similar percent level as in the case of the strong FSI) can be neglected when scaled by its 
contribution ~ 2r*/a. It concerns also the case of hadronic atoms since the subleading 
part is the same as in the continuous spectrum at k* = 0. 

Adopting the smoothness and equal-time approximations (with the accuracy of a 
few per mil), we can rewrite ( |69|) and (172|) for the production cross sections of particles 1 
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and 2 in continuous and discrete spectrum at low relative or binding energies as follows: 
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'd 3 P 6 v 7 ,1,z d 3 Pl d 3 p 2 - 

where 6 = {nO} and pi = P&mj/ (mi + m 2 ) in ( 1751) ; for equal-mass particles Pi = p 2 = 
P;,/2 and 7i = 7 2 = 7b- Particularly, for 7t + tt~ production, one can then rewrite the 
correction factors in ([38]) and ([39]) as 



l + 5(k*)^{\^v*)\%[A c (r])) 
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(77) 
(78) 



We will show that the r*-dependence of the wave functions f/'-k* an d V'f f° r two 
oppositely charged particles in continuous and discrete spectrum is practically the same 
at separations r* in the pair cm. system much smaller than the Bohr radius \a\. 
Therefore, the corrections to (175|) and f[76[) (arising due to the smoothness and equal-time 
approximations used in their respective derivation from ( 1121) and ( [20]) ) practically cancel 
out in the ratio of the numbers of pairs produced in continuous and discrete spectrum 
provided (r*) SL <C \a\. 



4-3.3. The effect of residual charge The formalism of section 2 assumes a free motion of 
a given particle pair during the final stage of the collision. Here we will estimate the FSI 
effect of the residual charge which is known to substantially influence particle spectra 
and, to a lesser extent, also particle correlations in low energy collisions involving nuclei 
[23] . Since, at high energies, this effect can be expected of minor importance, we will 
estimate only its upper limit. 

Generally, instead of the two-particle Bethe-Salpeter amplitude \I/W^(xi, x 2 ), the 
correlation function is determined by the amplitude \1/W^ a ^(xi, x 2 ) representing the 
solution of a complicated multi-body problem, taking into account interaction between 
the two particles and also their interaction with the residual system described by the 
quantum numbers {a}. For our purpose, it is sufficient to approximate these quantum 
numbers by an effective (comoving with a given pair) point-like residual charge Ze and 
consider a thermal motion of the two particles with the temperature T ~ m w in the rest 
frame of this charge. 

Let us start with the hypothetical case of particles that interact with the charge 
Ze but their mutual interaction is "switched off". In such a situation we can treat 
the systems (1,Z) and (2,Z) independently. Then the interaction with the Coulomb 
center just leads to the substitution of the spatial parts of the plane waves e ipiXi 
by the usual Coulomb wave functions: e~ ip4ri — > e~ iPir '$^. z where $^. z (rj) = 
e" 5 * ^A c (rji)F(-ir]i, l,ipi), p { = p^ + p^, ^ = (p^) -1 , a { = (uiZiZe 2 )' 1 is the Bohr 
radius of the system (i, Z) (taking into account the sign of the interaction) generalized 
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to the relativistic case by the substitution — > Ui of the particle masses by their 
energies, <5« is the Coulomb s-wave phase shift, A c (j]i) is the Coulomb penetration factor 
and F(a, l,z) is the confluent hypergeometrical function; see ( 1241) . ( l89i) and ( |90l) . For 
the complete amplitude we have: 

*£Jf(*i, *a) = e^^+^-^f ( ri )^f (r 2 ) 

Note that a small contribution of spin-dependent electro-magnetic forces is neglected 
here so that <f/( + ) sz = vp(+) z is independent of the total spin S of the particle pair. 

Let us now "switch on" the interaction between particles 1 and 2. Since we consider 
the relative motion of the two particles at characteristic distances much slower compared 
with their motion with respect to the Coulomb center, it is natural to assume that in 
such a case the plane wave e~ lk * r * in (179|) will be basically substituted by the Bethe- 
Salpeter amplitude ip~(x) describing the relative motion of isolated interacting particles. 
After this substitution we get the amplitude in so called adiabatic (factorization) 



approximation 

= e iPX ^f(x)$- z (rO$ p f (r 2 ). (80) 

Instead of the six-dimensional correlation function 7£(pi,p 2 ) we calculate the one- 
dimensional one, 7Z z (k*), with the numerator and denominator integrated over the 
simulated particle spectra. In the equal-time approximation, 

N(k*) 

£ E5P5|^ k .(r*)$^(r 1 J$ p f(r 2 J| 2 

**(*•) = ^ Nm ' : — , (si) 

£ |^(riJ^f(r 2i )P 
i=i 

where N(k*) is the number of generated particle pairs in a given k* bin. To separate 
the pure effect of the residual Coulomb field on particle correlations, we compare the 
correlation function TZ z {k*) with the one, TZ" Z '' (k*), taking into account for the latter the 
effect of the nucleus Coulomb field on one-particle spectra but not on particle correlations 
(i.e., simulating the argument r* independently of the arguments r x and r 2 ). Note that 
due to the velocity dependence of the correlation function, 71 " z " = 7Z Z=0 only at a fixed 
pair velocity v. In figure we present the ratios of the 7r + 7r~ correlation functions 1Z Z 
and 1Z z " assuming that the pions are emitted in the rest frame of the residual charge Z 
according to the thermal law with a temperature of 140 MeV at the space-time points 
distributed according to a product of Gauss functions with the equal spatial and time 
width parameters ro = ctq. One may see that even for the radius ro as low as 2 fm 
the effect of the residual comoving charge as large as Z = 60 is less than a few per 
mil. Taking into account that the effective radius ro is larger than 2 fm even for proton 
collisions with low-Z nuclei and that the effective residual charge is only a fraction of 
the target nucleus charge, one can conclude that the effect of the residual charge is on 
a negligible level of a fraction of per mil. 
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4-3.4- Femtoscopy with non-identical particles The FSI effect allows one to access the 
space-time characteristics of particle production also with the help of correlations of non- 
identical particles. One should be however careful when analyzing these correlations in 
terms of simple models like those assuming the Gaussian space-time parametrization of 
the source. The simplified description of the r*-separations can lead to inconsistencies 
in the treatment of QS and FSI effects. While the QS and strong FSI effects are 
influenced by large r*-separations mainly through the correlation strength parameter 
A, the shape of the Coulomb FSI is sensitive to the distances as large as the pair Bohr 
radius (hundreds of fm for the pairs containing pions). 

This problem can be at least partially overcome with the help of imaging techniques 
[4"0] or transport simulations. The former yield the redistribution inverting the 
measured correlation function using the integral equation (j75l) with the kernel given 
by the wave function squared. The latter account for the dynamical evolution of the 
emission process and provide the phase space information required to calculate the QS 
and FSI effects on the correlation function. 

Thus, the transport RQMD v.2.3 code was used in a preliminary analysis of the 
NA49 7r + 7r~, 7i + p and ix~p correlation data from central Pb+Pb 158 A GeV collisions 
[20J. The model correlation functions 7£rqmd(Q; s r ) have been calculated using the FSI 
code based on the formalism developed in [12] , weighting the simulated pairs by squares 
of the corresponding wave functions. The scale parameter s r , multiplying the simulated 
space-time coordinates of the emitters, was introduced in the model correlation function 
to account for a possible mismatch of the redistribution. For this, a set of correlation 
functions 71rqmd(Q] s l r ) was calculated at three chosen values s r of the scale parameter 
and the quadratic interpolation was used to calculate 7?.rqmd(<3; s r ) for arbitrary value 
of s r : 

^•rqmd(Q; s r ) = 77 J777 ^-^-rqmd(Q; si), (82) 

t=l [ s r ~ Sr )\ s r ~ s r) 

where {i,j,k} are permutations of the sequence {1,2,3}. The NA49 correlation 
functions were then fitted by 

K{Q) = N [\TZ RQMD (Q; s r ) + (1 - A)] (83) 

with two additional parameters, the normalization N and the correlation strength A. 
The fitted values of the A-parameter are in reasonable agreement with the expected 
contamination of ~ 15% from strange particle decays and particle misidentification. The 
fitted values of the scale parameter show that the RQMD transport model overestimates 
the r*-separations of the pion and proton emitters by 10-20% thus indicating an 
underestimation of the collective flow in this model. 

The shape of the correlation function is less influenced by large r*-separations in the 
case of two-particle systems with the absent Coulomb FSI, e.g. in the case of pA system. 
The data on pA correlations in heavy ion collisions show a significant enhancement at low 
relative momentum, consistent with the known singlet and triplet pA s-wave scattering 
lengths. In fact, the fits using the analytical expression for the correlation function [12] 
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yield the Gaussian correlation radii of 3-4 fm in agreement with the radii obtained from 
pp correlations in the same experiments. These radii are smaller than those obtained 
from two-pion and two-kaon correlation functions at the same transverse momenta [H] 
and are in qualitative agreement with the approximate m t scaling expected in the case 
of the collective expansion, see (1671) and ( 1681) . 

The correlation function of non-identical particles, compared with the identical 
ones, contains a principally new piece of information on the relative space-time 
asymmetries in particle emission [12]. Since this information enters in the two-particle 
FSI amplitude through the terms odd in k*r* = p*(r* — r|), it can be accessed studying 
the correlation functions lZ + i and IZ-i with positive and negative projection k* on a 
given direction i or, - the ratio 1Z + i/1Z_i. For example, i can be the direction of the 
pair velocity or, any of the out (x), side (y), longitudinal (z) directions. In LCMS, we 
have r* = r iy except for r* = Ax* = 7*(Ax — v t At), where jt and v t are the pair LCMS 
Lorentz factor and velocity. One may see that the asymmetry in the out (x) direction 
depends on both space and time asymmetries (Ax) SL and (At) SL . In case of a dominant 
Coulomb FSI, the intercept of the correlation function ratio is directly related with the 
asymmetry (r*) SL scaled by the pair Bohr radius a: 



It appears that the out correlation asymmetries between pions, kaons and protons 
observed in heavy ion collisions at CERN and BNL are in agreement with practically 
charge independent meson production and, assuming m\ < m 2 , with a negative 
(Ax) SL = (x± — £2) SL and/or positive c(At) SL = c{t\ — t2) SL on the level of several 
fm [201113]. In fact they are in quantitative agreement with the RQMD transport model 
as well as with the hydro-motivated blast wave parametrization, both predicting the 
dominance of the spatial part of the asymmetries generated by large transverse flows. 

In the thermal approach, the mean thermal velocity is smaller for heavier particle 
and thus washes out the positive spatial shift due to the flow to a lesser extent. As 
a result, (x n ) Sh < (xx) SL < (%>) SL - The observation of the correlation asymmetries 
in agreement with the mass hierarchy of the shifts in the out direction may thus be 
considered as one of the most direct signals of a universal transversal collective flow 
[20] . This is in contrast with the effect of m t scaling of the correlation radii which can 
be also explained by a large transverse temperature gradient like in the Buda-Lund 
model [H]. 

4-3.5. Correlation measurement of strong interaction One can also use the correlation 
measurements to improve knowledge of the strong interaction for various two-particle 
systems. In the collisions involving sufficiently heavy nuclei, the effective radius r 
of the emission region can be considered much larger than the range of the strong 
interaction potential. The FSI contribution is then independent of the actual potential 
form [15]. At small Q = 2k* and a given total spin S, it is determined by the s-wave 
scattering amplitude f s {k*) [12]. In case of \f s \ > ro, this contribution is of the order 
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of \f s /r \ 2 and dominates over the effect of QS. In the opposite case, the sensitivity of 
the correlation function to the scattering amplitude is determined by the linear term 
f S /r . 

The possibility of the correlation measurement of the scattering amplitudes has been 
demonstrated [20] in a preliminary analysis of the NA49 7r + 7r~ correlation data within 
the RQMD transport model. For this, besides the r*-scale s r , the strong interaction 
scale Sf has been introduced in the RQMD correlation function 1Z(Q; s r ,Sf), rescaling 
the original s-wave 7r + 7r~ scattering amplitude taken from |5j: f(k*) — > Sff(k*); it 
approximately corresponds to the rescaling of the original scattering length f Q = 0.232 
fm. The fitted parameter Sf — 0.63 ± 0.08 appears to be significantly lower than unity. 
To a similar but somewhat weaker rescaling (~ 0.8) point also the preliminary result of 
experiment DIRAC on the pionium lifetime [4] , the BNL and CERN data on K\± (6] and 

— > n^^n [7] decays as well as the two-loop calculation in the chiral perturbation 
theory with a standard value of the quark condensate [8]. 

Comparing the fit results with the theoretical predictions, one should have in mind 
that the latter are subject to the electro-magnetic corrections on the level of several 
percent and that the correlation measurement may underestimate f(k*) by a few percent 
due to the use of the equal-time approximation. A substantial systematic error can also 
arise from a simplified fit of the strong FSI amplitude. To avoid the latter, one can use 
the Roy equations and represent the 7r + 7r~ strong interaction amplitude at low energies 
as a unique function of the isoscalar and isotensor s-wave scattering lengths a[j and a^, 
see Appendix D in |46j. The two-parameter dependence of the scattering amplitude can 
be further reduced to a single-parameter one within the generalized chiral perturbation 
theory predicting a strong correlation between the two s-wave scattering lengths (see 
equation (13.2) in [8J). The systematic error due to the uncertainty in the fitted re- 
distribution (e.g., in the scale parameter s r ) can be diminished in a simultaneous analysis 
of 7r + 7r~ and 7r ± 7r ± correlation functions. The high statistics DIRAC data on two-pion 
correlations may thus allow one to determine the s-wave scattering lengths afj and 
better than to 10% and serve as complementary to the pionium lifetime measurement 
in the same experiment. 

The correlation technique was also used to estimate the singlet AA s-wave scattering 
length based on the fits of the AA correlation data from Pb+Pb collisions at 158 A GeV 
|20j. Though the fit results are not very restrictive, they likely exclude the possibility 
of a large singlet scattering length comparable to that of ~ 20 fm for the two-nucleon 
system. Similarly, the fit of the pA and pA correlation functions measured in Au+Au 
collisions at v^nn = 200 GeV allowed one to determine the corresponding spin-averaged 
s-wave scattering length. The fitted imaginary part of the scattering length of ~ 1 fm is 
in agreement with the pp results (thus pointing to about the same pA and pp annihilation 
cross sections) while the real part appears to be more negative [471] . 
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5. One-channel wave functions 

5.1. Continuous spectrum 

5.1.1. Short-range FSI Let us start with the case when the two-particle FSI is due to 
the short-range forces only. In the considered region of small k* the short-range particle 
interaction is dominated by s-waves. Since the radius of the s-wave interaction is usually 
small compared with the distance r* between the production points of particles 1 and 2 
in their cm. system, the FSI effect is mainly determined by the asymptotic behaviour 
of the scattered wave outside the region of the strong interaction r* > d: 

A^_ k ,(r*) = /(F)e ifcV 7r*. (85) 

The s-wave amplitude / depends on the magnitude of the vector k* only. Assuming the 
absence of inelastic transitions, it satisfies the one-channel s-wave unitarity condition 
= k*\f\ 2 or, equivalently = —k*, and so can be represented as 

/ = [exp(2i5 ) - l]/(2iF) = (K- 1 - iF)" 1 , (86) 

where S is the s-wave phase shift and K~ x = k* cotS is a real function of k*. Usually 
(for potentials vanishing with the distance exponentially or faster), this function is real 
also for negative kinetic energies k* 2 /(2fi), so that its expansion can contain only even 
powers of k* [25]. Retaining near threshold only the first two terms in the expansion, 
one can express the function K^ 1 or K through the corresponding two parameters: 
scattering length f and effective range d or curvature b : 

K- l = f l + ±d k* 2 K = f + b k* 2 b = -U f 2 . (87) 

The expansion of K~ x is superior for two-nucleon systems (due to large scattering 
lengths, amounting to about 20 fm in the singlet case) while for other systems, the 
.fT-expansion is often preferred. To extend the latter to a wider energy range, it is 
usually written in a relativistic form and additional parameters are added. For example 

K = ^ fth^o £ A 2j x = 2k y^ (g8) 
Vs s-s ^ 

where s = (pi + p 2 ) 2 = (u* + ul) 2 , cu* >2 = {m\ 2 + k* 2 ) 1/2 and s th = (m 1 + m 2 ) 2 . 
The parameter so takes into account the eventual resonance, specifying the value of 
the two-particle invariant mass squared where the phase 5o{k*) passes through 90°. The 
behaviour of the s-wave i^-function in a wide fc*-interval is however of minor importance 
since we are interested in the near-threshold region and, have already neglected the 
p-wave correction 0(k* 2 ai/r*) in (I8"5"j) ; here a\ is a p-wave scattering length. For 
7T + 7r _ system, a\ ~ 0.1 fm 3 , /o m 0.2 fm, do m —10 fm and the relative p-wave 
contribution to the k* 2 term due to the short-range FSI composes in the production 
cross section ~ cli/(cl\ — d f Q 2 /2 — /o 3 /3) ~ 30%; the relative contribution of the k* 2 
term ~ (a x — d f 2 /2 — /o 3 /3)/c* 2 //o being less than 1% of the total short-range FSI 
contribution for Q = 2k* < 30 MeV/c. 
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Table 1. The pair Bohr radius including the sign of the interaction, a = (/izi^e 2 ) , 
and the characteristic width of the Coulomb correlation effect, Q c = 2k* = 4-7r/|a|, 



correspondin; 


I to \ri\- 1 


= 2tt (see 


(pM)) and the first panel in 


figure O. 
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Note that the extension of the asymptotic wave function in the inner region leads 
to a relative shift in the production cross section of the order \f\ 2 jp^-Re(l//)/((r*) SL ) 3 
[12] HE]. The leading part of this shift can be, in principle, corrected for (see subsection 
6.3). However, being quadratic in the amplitude /, it is rather small for ittt-, ttK- or 
7rp-systems - usually not exceeding several percent of the short-range FSI contribution. 

5.1.2. Account of the Coulomb FSI Similar to the case of neutral particles, we will 
approximate (with the same accuracy) the wave function of two charged particles near 
threshold by the asymptotic solution outside the region of the strong interaction r* > d. 
It is well known that the long-range Coulomb interaction modifies both the plane and 
spherical waves [25] : 

^(r*)=e i ^ V / A^) 

where £ = k*r* + k*r* = p(l + cos6'*), p = k*r*, r\ = (k*a)~ l , a = (/xzi^e 2 ) -1 is the two- 
particle Bohr radius including the sign of the interaction (see tabled]), 5 C = argr(l + ir)) 
is the Coulomb s-wave phase shift, A c {rj) is the Coulomb penetration factor defined in 

ii, 

F(a, l,z) = l + az/V. 2 + a(a + l)z 2 /2\ 2 + ■■■ (90) 

is the confluent hypergeometric function and G = -JA~ c (Gq + iF ) is a combination of 
the regular (F ) and singular (G ) s-wave Coulomb functions (see, e.g., [4"5]): 

G(p, 77) = P(p, ri) + 2 V p [In \2 V p\ +2C-1 + X (v)] B(p, rj). (91) 

Here C = 0.5772 is the Euler constant, the functions 

00 

B(p,r])=J2Bs B = l B 1 = V p 

s=0 

00 

P(p, V )=J2Ps Po = l Pi = Q ••• (92) 

are given by the following recurrence relations: 

(n + l)(n + 2)B n+1 = 27] P B n - p 2 B n _ 1 

n(n + l)P n+1 = 2 V pP n - p 2 P n ^ - (2n + l)2 V pB n , (93) 



-ik*r* 



F(-iT],l,iZ) + f c (k* 



(89) 
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B(p,rj) = F /(/>v/^Q ~~ *■ S i n (p)/P an d P(p,v) ~> cos(p) in the limit 7/p = r*/a — > 0. 
The function 

X (r ? ) = /i(r/) + iA c (r ? )/(2r ? ), (94) 

where the function h(r]) is expressed through the digamma function ip{z) = T'(z) /T(z): 

h{rj) = U>(irj) + ip(-iv) - ln(r? 2 )] A (95) 

For \r]\ < 0.3 the function h{rj) = \.2if — In \r\ — C, while at large \r]\ this function can be 
represented by a truncated series in inverse powers of rj 2 : h(r]) = ?7 — 2 /12 + 77" 4 /120 + - ■ •. 
The amplitude 

fc(k*) = f(k*)/A c ( V ), (96) 

where f{k*) is the amplitude of the low-energy s-wave elastic scattering due to the short- 
range interaction renormalized by the long-range Coulomb forces. Assuming again the 
absence of inelastic transitions, the amplitude f{k*) = (e 2l<5 ° — l)/(2iA;*) and satisfies 
the one-channel s-wave unitarity condition. Similar to the case of neutral particles, one 
then has [25] : 



/c(**) = (k- 1 - 2 -^\ 1 , (97) 

where the function K can be parametrized according to ( J87l) or (!88l) . 

Note that 5 C -»• 0, A c -> 1 for r, -> (k* > \cl\- 1 ) and G -> e^, F -> 1 for 
r/p = r*/a — > 0. So, the two-particle wave function in the absence of the long-range 
Coulomb forces is recovered provided r*, /o and are much smaller than the Bohr 
radius \a\. 

In figure El we plot A c (r)) and x(^) as functions of the variable \ri\~ 1 = \ak*\. For 
the system of two charged pions, this variable approximately corresponds to Q — 2k* 
in MeV/c. At k* — > 0, the Coulomb penetration factor A c (rj) respectively tends to 
and 00 for like and unlike particle charges. With the increasing k*, this factor slowly 
approaches unity: A c (i]) pa 1 — ttt] for k* > 2ir/\a\. Note that the quadratic behaviour of 
^x(v) = h{ff) pa i]~ 2 /l2 at |?7| _1 < 1 is changed by a steep quasi-linear rise in the interval 
1 < l^l -1 < 5; the corresponding slope being about 0.26. As for ^sx(v) = A c {rj) / \2rj) , at 
k* = it equals to and —it for like and unlike charges respectively, and, for k* > 27r/|a|, 
it approaches the linear /^-dependence: ^x(v) ~ (v 1 ~ 7r )/2- 

5.1.3. The small- and large-r* limits Since we are interested in the region of small 
relative distances r* compared with the Bohr radius \a\ and small relative momenta 
Q = 2k* compared with 1/r*, it is useful to write the first terms in the expansion of the 
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hypergeometric functions F and G in r*/a and p 
F(-i7?,l,i£) = ! + -(! + x) 



k*r*. We have (x = cos 9* 



l + ^n+x)- 



p 



;i + x) 2 + o( P 3 



+ o 



,r 
a 



G(p, V ) = l 



2 a 



(98) 



In |2 — | + 2C- l + x(v) 
a 




2 ^ r *2 one can 



For some systems of interest (tttc, ttK, up), \f Q \ 2 
neglect the Q-dependence of the scattering amplitude and, after the averaging over the 
uniform x-distribution, write the correlation function at a fixed separation r* as 

(l^-k*(r*)| 2 ) 

Jo 



K(k*;r*] 
= Mv) 




ik*r* 



F*G 



Ml) 



r* f f 
f 2— + 2— + 2— 

a r* a 

9 a 3 r* I 




1 + h{rj) 



(99) 



+ 0(p 4 



Note that in the case of an anisotropic r*-distribution, Eq. (199]) implies the integration 
over the direction of the vector k* = Q/2, distributed isotropically for non-correlated 
particles at Q — > 0. In the case of the cut Qt < Q™* on the component of the vector Q 
transverse to the direction of the pair three- velocity v and Q > Q^, Eq. (I9"9l should 
be modified by the substitution 2/9 — ► g cut 2/9, 



Pcut - 1 — 2 (|(3cxxs y ^ - 1) rf )jiu ! 6 (0.0.2 



(100) 

(Qr/Q) 2 ] 1/2 



IS 



where 6 r * is the angle between the vectors r* and v, and c m - m = [1 
the minimal absolute value of the cosine of the angle between the vectors k* and v. For 
pion pairs containing pions from resonance decays, one may expect (r£) > (r^) |49j (i.e. 
(cos 2 # r *) > 1/3) and so g cut > 1. 

In figures El M and M, we show the Q-dependence of the functions B(p,rj), 
P(p, rj), G(p, if) and the reduced correlation function 1Z/ A c as well as the corresponding 
main contributions due to the interference term and the modulus squared of the 
hypergeometric function for the tt + tt~ system at r* = 5, 15, 50 fm. One may see 
that the almost universal quasilinear decrease of TZ/A C for r* <~ 20 fm is due to the 
interference term, and that it is changed, for higher r*-values, by a steep rise due to 
the |F| 2 -term. It appears that the linear fit of 7Z/A C recovers the intercept better than 
to 2 per mil for r* <~ 20 fm and - better than to 2 percent even for r* = 50 fm (see 
table E]). 
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Table 2. Results of the linear fits of the reduced ir + ir~ correlation function: 
1Z/A C = cq + ciQ in different intervals < Q < Q m ax- The function 1Z/A C is calculated 
at r* — 5, 15 and 50 fm in the approximation of a constant scattering amplitude 
fc(k*) = /o = 0.232 fm and, assuming the uniform distribution of the cosine of the 
angle between the vectors r* and k* = Q/2. Also shown are the corresponding values 
of TZ/A C at Q = (the intercepts). 



r* , fm 


Intercept 


Q max , MeV/c 


10 


20 


30 


40 


50 


5 


1.077 


Cf) 




1.077 


1.077 


1.077 


1.077 


1.078 






Cl, 


(GeV/c)- 1 


-0.55 


-0.47 


-0.48 


-0.52 


-0.57 


15 


0.961 


CO 




0.961 


0.961 


0.961 


0.960 


0.959 






Cl, 


(GeV/c)- 1 


-0.59 


-0.56 


-0.55 


-0.51 


-0.42 


50 


0.783 


C{) 




0.778 


0.768 


0.766 


0.773 


0.783 






Cl, 


(GeV/c)- 1 


2.55 


4.61 


4.99 


4.38 


3.69 



To clarify the origin of the quasilinear behaviour of the reduced correlation function 
1Z/A C , one can use ( |99l) to estimate the slope at small Q: 

fH\'_ d (K\ ., of dh Mr* 2f \ , 2 
UJ = dQ [aJ = ±2/ W " 1 9 R + 3 V*) T Q ' (101) 
where the sign + (— ) corresponds to the Coulomb repulsion (attraction). Using the fact 
that dh/d\ri\- 1 « 0.26 for 1 < < 5, one has (K/A c )' « -(0.6 + bQ) (GeV/c)" 1 for 
the 7T + 7r _ -system at 1 < Q < 5 MeV/c and fo = 0.232 fm, where b is small (b <~ 0.03 
(MeV/c)" 1 ) and positive for r* <~ 20 fm and, for larger r*-values, b is negative and 
its magnitude rapidly increases with r*. As a result, the slope of the reduced 7r + 7r~ 
correlation function is negative in this Q-interval and nearly constant for small re- 
values, while it becomes positive and rapidly increases with Q for r*-values of several 
tens of fm or larger. For Q > 5 MeV/c, the absolute value of the slope due to the 
/i-function decreases as ~ 2.35|r/|. It appears that, for the 7r + 7r~-system at r* <~ 20 
fm, this decrease is approximately compensated by the Q-dependence of the functions 
B, P and F, (i.e., at p < 1, by the second term in (jlUIjl ). so that (R/A c )' ta -0.5 
(GeV/c)" 1 up to Q = 50 MeV/c. 

Note that the Qr-cut substantially influences the Q-dependence of the reduced 
correlation function only for sufficiently large values of r*/\a\, leading to the substitution 
r*/|a| — > ^cut^*/|a| in (HO 11) at Q > Qt ■ Particularly, for 7r + 7r~ pairs containing an uo- 
decay pion, one may expect g cut > 1 and so a more steep rise of the reduced correlation 
function at Q > Qt- 

To estimate the behaviour of the correlation function at large r* or k*, one can 
exploit the known asymptotic expressions for hypergeometric functions. Thus, at 
f > 1 +r/ 2 , 

v / A c (r/)F(-ir/, 1, if) -> M - i^ e^+^t) + !L e Ks c +^<) (1Q2) 
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and, at p ^> 1 + rf, 

G(p,rj) - ^[AJ^)e^ + P-^\ (103) 

so that both the effects of the Coulomb and strong FSI vanish in the cross section as 
r*~ 2 . In fact, the asymptotic expression for the F-function in ( 11021) cannot be used in 
the case of nearly opposite directions of the vectors k* and r* (cos#* ~ —1) when the 
variable £ = p(l+cos#*) is suppressed even at large p = k*r* . This leads, after averaging 
over the angles, to a slower vanishing of the Coulomb effect, as r* -1 , in agreement with 
the classical Jacobian factor [1 - 2/ (ar*k* 2 )] 1/2 m 1 — (ar*k* 2 )~ l . 



5.2. Discrete spectrum 



5.2.1. General s-wave solution Since the Schrodinger equation at a small negative 
energy — e& = — k 2 /{2p) practically coincides with that in continuous spectrum at zero 
energy, the r*-dependence of the corresponding wave functions at given orbital angular 
momentum I and r* <C k^ 1 is the same. This important conclusion was first stated by 
Migdal for the pn-system [26]. In fact, both solutions (at positive and negative energies) 
can be written in the same form for any r*, up to an energy dependent normalization 
factor N '. Outside the region of the short-range interaction, r* > d, we can write the 
s-wave solution as a combination of the regular and singular Coulomb functions: 



^=o(r*)=A%) 

.Py/Mv) 

At d < r* <C \a\ and \p\ <C 1 it takes on the form: 



F o{p,v) | f (k*) G ( p,V ^ 



(104) 



iMO = + -) + o((-) 2 ) + o( P 2 ) + 

a \ a / 



fc 



1 + 2— [ In 12— I + 2C 



a 



X (l + -) + (- 



0(p 2 



(105) 



For positive energies, M = e ldc yjA c (i]) and, at k* — > 0, f c = f /[l - 2f x(±oo)/a] 
x(+oo) = (a > 0) or x(— oo) = — in (a < 0). In the case of opposite charges (a < 0), 
( TT05D yields: 



1=0 



r ' 
\a\ 



fo 



In 



2r* 



+ 2C 



fo 



+ 



V r* , 



2ivrA + 0((-) 2 )+0(p 2 ) 



(106) 



For the discrete levels at negative energies, the substitution k* 
particularly yielding [25l 148] : 

1 N 



7T 7T 1 

x( ^ ) = 2 cot( ^y ) + 2 



21n(ft n |a|)+V 



m n has to be done, 
1 



/'b^ | | 



(107) 



where 



I 1 fa n Qj 



-1 



A more compact form of fjlOTj) follows from the relation 



t/)(— x) = if)(x) + 7r cot(7ra;) + x l : 



x{Vn) = vrcot(7rx„) - (2x n ) x [0(x„) - 3] x n = (n n \a\) 1 



108) 
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(x) = 2 + 2x[lnx — ip{x)]. 



(109) 



5.2.2. Energy levels For a pure Coulombic atom (a < 0, / = 0), only the solution 
F /p, regular at r* — > 0, contributes and the requirement of its exponential damping at 
large distances fixes the energy levels. The corresponding /t-values at a given principle 
quantum number n are equal to K c n = (n|o|) -1 . The wave functions ip, 



coul / 
nl \ 



can then 



-22+1 



be expressed in terms of the Laguerre polynomials L n+l _ lK 

' 2r* 



For / = 0, 



coul / 
nO \ 



7 /,coul 



(0)exp 



Li 



n-l > 



[n- nl. 



(110) 



n\a\j \n\a\ 

The square of the wave function ip^ l (0) at zero separation is given in ff25l) and the 
Laguerre polynomials are defined by the following recurrence relations: 

n-l 

£i_!(*) = (n.n!)5:d(*) 

(111) 



s=0 



At r* -C nlal 



coul / 
nO l 



n°0 Ul (0) 



s(s + 1) 



i-n + o (-) : 



(112) 



The strong interaction slightly shifts the Coulombic energy levels thus making the 
regular part of the general solution (I104j) divergent at large distances. Therefore, the 
amplitude f c has to have a pole at k* = iK n , and so, according to ( 1971) . 

fodo 



XiVr, 



l + O 



2K(tK n ) 2/ 

Using (11081) and (11131) , one can fix the energy levels E, 
with the relative error of 0(a -3 ): 



na) 



(113) 

K n 2 / (2/i) in discrete spectrum 




1 + /o< 



n 



i] 

n|a|) _1 . 



4tt z 



(114) 



To show this, one can put n n = 
equality tan(7nr n ) = — tan(7rx„e) 
\d)(x n 



X r 



n/(l + e) and use the 



<(l + e) 

= — (7ra; n e)[l + (7r 2 n 2 /3)0(e 2 )] and the inequality 
4>{n)\ < O(e), the latter following from the fact that (f)'(n) vanishes faster 
than n~ l . Equation (11141) is in agreement with the result of [50J for the relative energy 
shift e(n, 0) = (2 + e)e = e (n, 0)[1 + e (n, 0)pi(n, 0)], where e (n, 0) = 4/ ^ and 
Pi(n, 0) = 0(n)/4. The function 0(n) is defined in ( I109f) with the digamma function for 
the integer values of the argument given by the recurrence relation: 

if>(n + l)=ip(n) + l/n = -C = -0.5772. (115) 

Note that </>(n) w 3 is nearly constant: 0(1) = 2 + 2C = 3.15443, 0(2) = 3.08145, 
0(3) = 3.05497, . . ., 0(oo) = 3. 
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5.2.3. Normalization Since N(j] n ) = (to compensate for the pole of the amplitude 
f c at k* = iK n ), the s-wave solutions in discrete spectrum are now given (for r* > d) by 
the second term in (11041) . exponentially vanishing at large distances: 

G(p n ,r] n ) 



M'{n)K{iK n ) 

G(p n ,T) n ) 



1 + 01 



fodp 
n 2 a 2 



(116) 



The arguments p n and r] n are taken at k* = in n and the normalization factor 

N'{n) = M{ri n )f c {iK n )/K{iK n ) 
is set by the requirement 

|io(rfd¥ = l. 



;ii7) 



(118) 

Note that the extension in the integral fll 18j) of the asymptotic wave function (11161) into 
the inner region r* < d leads to negligible relative errors O(f d 2 /(net) 3 ), 0(/qC?/ (na) 3 ) 
in the normalization factor TV'. Using the expansion of the G-function in the square 
brackets in (11051) and the expression for xiVn) i* 1 <1H3fl . one can write for distances 
d < r* < \a\: 




+ 2C-1 



fa 



;il9) 



Comparing ( 11191) with the low-r* expansion ( j 1 121) of the pure Coulombic wave function 
and, also taking into account the exponential damping at large distances, one can 
approximate the wave function (11161) at r* <C \a 2 /fo\ by the expression: 

N'(ri) 



4>. 



app/ 
nO I 



^°o ul (0) 



coul / 
nO ' 



1-2A 


fan 


2r* 






\a\ 




a 







(120) 



From the results of calculations for the s-wave 7r + 7r~ atoms, presented in upper panel 
of figure [TOl one can see that the squares of the approximate and exact expressions 
(11201) and (11161) practically coincide for the distances up to several tens fm and that the 
agreement is better than percent even at r* ~ \a\. 

It follows from (11201) that the relative difference of the normalization factors N'(ri) 
and ^o ul (0) scales as O(fo/ a )- In fact, this difference can be fixed when extending the 
theory to a multichannel case and requiring the equality of the total width r n = —2QE n 
and the sum of the partial widths (see (1111!) and (fR5l) or, (!R7l) . fjR8|) and (jRij) ). As 
a result: 

2/o 



|A^(n)/Co Ul (0)| 



1 



n 



121 



n\a\ 



We have checked f 1 1 2 1 [) . calculating M' from the integral (11181) for various values of the 
scattering length f , Bohr radius \a\ and the principle quantum number n. 
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5.3. Universality 



Comparing ( 11061) and ( 11191) . valid for the distances d < r* <C \a\, one confirms 
the important conclusion, already stated at the beginning of subsection 5.2, about 
the universality of the r*-behaviour of the moduli squared of the s-wave solutions in 
continuous (k* — > 0) and discrete spectrum, up to corrections vanishing as inverse 
squares of the Bohr radius \a\. Assuming f <~ d, one has: 



A 



n0\ 



^ (r*)/«(0) 



- 1 



f 2< 



01 




(122) 



0(P 2 



The universality holds with the same accuracy also if the s-wave solution in continuous 
spectrum were substituted by the complete wave function (recall that V-°k*(0) = 
ip c k °4{0) = Al' 2 ), provided the averaging over the angle between the vectors r* and 
k*: 



V>-k*(r*)/«(0) 



^ n o(r*)/Af'(n) 



1 



4tt 2 OI 



7o 



o 



fpdp 
n 2 a? 



,*2" 



(123) 



+ 01 



+ 0( P 2 



This result follows from the fact that, at k* — > and typical distances r* <C \a\, the 
total wave function in continuous spectrum almost coincides with the s-wave amplitude 
V'A:*o( r *) ( see the lower panel in figure [9]): 



?A_ k *(r*) =4> k * (r*) +e 



k*r* 

k*a 



O - +0(p 2 



(124) 



and, that the relatively significant correction term 0(r*/a) in the square of the wave 
function ?/Ak*(r*) vanishes after averaging over the direction of the relative three- 
momentum Q = 2k* or, - after suppressing the signs of the components Qi of the 
vector Q (assuming a symmetric detector acceptance with respect to the reflection 
Qi —> —Qi)- From the lower panel of figure [TUl one can see that for the 7r + 7r~ system, 
the universality holds to better than percent for r* <~ 50 fm. Note that A° (not 
shown in figure [TUl) is negative and, contrary to A°, it shows the strongest deviation 
from zero for n = 1, achieving a per mil level already at r* w 20 fm. 

Comparing (!75|) and (!76|) . one can see that the number Na of produced 7r + 7r~ atoms 
is determined by the number of non-atomic 7r + 7r~ pairs in the region of small k* . So, Na 
is actually proportional to the ratio (|V ; no( r *)| 2 )op/(|V'-k*( r *)| 2 )gp in which the effects 
of the r*-dependence as well as the corrections due to non-equal emission times (t* ^ 0) 
and smoothness assumption are to a large extent compensated for, being practically the 
same for the wave functions in continuous spectrum at k* ^ and discrete spectrum 
at r* <C \a\. In fact, according to ( 1121111231) . one can write the ratio of the finite-size 
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correction factors at small relative momenta (k* <C (l/r*) SL ) and moderate distances 
between the particle emitters ((r*) SL <C \a\) as 



1+jn 
1 + 



(l^no(r*)/^°o ul (0)| 2 l^ 

2\SL 



(lfk.(r*)/CS' 



1 + <£(n 



2/o 



n\a\ 



O 



/ r *2\Sl/ 



iV(n) 



Co ul (0) 



1 + <5(fc*) 
+ 0(F 2 (r* 2 ) SL )| 



(125) 



7T 1 7T 



sec 



thus leading to (|42p up to a small correction due to the transition 7r°7r° 
(11541) ). Recall that though the /c*-dependence of the correction factor in braces is 
quadratic at very low values of k*, in fact, in a wider k* -interval and for sufficiently 
small values (r*) SL <~ 10 fm, it shows a quasi-linear and almost universal behaviour 
(see figure [9] and pQ ) . 



5.4- The n-dependence 

Neglecting the production of the tt + tt~ atoms with the orbital angular momentum I > 0, 
suppressed by powers of the ii + 7i~ Bohr radius \a\, the pionium production probability 
is given in (I34"j) and depends on the main quantum number n as 

w n oc (1 +5 n M°o ul (0)| 2 °c (1 +Sn)/n\ (126) 
The correction factor (1 + 5 n ) slightly modifies the n -3 law of simple ansatz in (123]) . It 
follows from ( 11251) that the n-dependence of the short- distance part of the correction 
S n is dominated by the renormalization effect of the strong FSI on the two-pion atomic 
wave function and that the renormalized correction (see also (142]) ) 

5' n = 5(0) + O{(r* 2 /a 2 ) SL ). (127) 

The most right panel in figure [12] confirms that the short- distance part of the correction 
S' n is practically independent of n and equal to 5(0). The renormalization correction 
2cj)(n)f /(n\a\) ~ 6f /(n\a\) is the largest for low values of n. For example, for pionium 
at n = 1 it composes ~ 0.3% (see also the central panel of figure [10]). As for the uj 
and rj' contributions to 8' n , their n-dependence is not negligible and the shifts from 5(0) 
compose up to ~ —0.004 and ~ —0.1, respectively. 

In [5T1 152]. the effect of the strong interaction on the n-dependence of the pionium 
wave function has been studied numerically, solving the corresponding Schrodinger 
equations. Thus, in pjT] . the ratio R n = ipno/ip^ 11 and the difference AR n — R\ — R n 
have been calculated for n = 1 — 3 using an exponential form of the short-range potential. 
According to (11121) . (11191) and (11211) . one has, up to corrections O(/ /a) and 0(r* 2 /a 2 ): 




From figure 1 of [51], one can deduce a value of ~ 0.15 fm for the scattering length f to 
achieve an agreement with the prediction of (11281) for the ratio R n at d < r* \a\. The 
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differences AR n , presented in figure 1 of [51 J for n = 2 and 3, are however by a factor 
of 1.6 higher than the corresponding predictions of (j!28p . For example, for 10 3 Ai? n at 
r* = 8 fm, n = 2 and 3, one can read from this figure the values^ 1.0 and 1.3 while, (11281) 
respectively predicts 0.6 and 0.8. This discrepancy may indicate that the calculation 
error, declared in [51] to be better than 10~ 4 , was underestimated by a factor of 5. 

In [52], a more refined numerical study of the n-dependence has been done 
accounting for the second channel (7r°7r°) and extended charges. The hadronic titt 
potentials have been chosen to reproduce the phase shifts given by two-loop chiral 
perturbation theory. The quantity d n = n 3 / 2 ip n0 /ip 10 —l has been calculated for n = 1—4. 
Similar to ( 11281) . one has for d < r* <^ \a\ 



up to corrections O(for* /a 2 ) and 0(r* 2 /a 2 ). The results of numerical calculations 
presented in figure 2 of [52] are in qualitative agreement with ( 11291) . d n being almost 
constant (except for the region of very small r*) and showing the right n- dependence: 
d n ~ — (1 — 1/n). Similar to [51], the numerical results for \d n \ are however higher, now 
by a factor of 2.5, than the predictions of (I129[) calculated with f = 0.2 fm which should 
correspond within ~ 10% to the choice of the potentials in [52J. Since the presence of 
the second channel leads to a negligible modification of (11291) (/o — > ?R,A aa w /o; see 
next chapter) and the correction due to the extended charges is also expected to be 
negligible (~ — |(r 2 )„-/a 2 ), the discrepancy in the size of the correction d n has to be 
attributed to the insufficient calculation accuracy or, to the incorrect matching of the 
scattering length. 

6. Two-channel wave functions 

6.1. Continuous spectrum in both channels 

It was implied until now that a long-time FSI takes place and can be separated 
in the Bethe-Salpeter amplitudes in the near-threshold final-state elastic transitions 
1 + 2 — > 1 + 2 only. In principle, however, it can be separated also in the inelastic 
transitions, 1 + 2 — > 3 + 4, characterized by a slow relative motion in both entrance and 
exit channels. The necessary condition for such a separation is an approximate equality 
of the sums of particle masses in the intermediate (777.3 + W4) and final (mi +777.2) states. 
Some examples are the transitions 7r + 7r~ <-> 7r°7r°, n~K + <-> tt°K or n~p <-> 7T°n. 
For such processes only the second term in the upper diagram in figure [1] contributes, 
now with the particles 3,4 in the intermediate state. In the equal-time approximation, 
the corresponding amplitudes reduce to the wave functions describing a two-channel 
scattering of the particles 1, 2 with the inverse direction of the relative three-momentum: 

§ One should correct the figure by interchanging the curves. The author is grateful to 
O. Voskresenskaya for pointing out this misprint. 




(129) 
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Table 3. The /3-channel momenta kp calculated at the a-channel thresholds 
fc* = 0. Also shown are the relative shifts Akp/kp arising from the non-relativistic 
approximation in the second formula in (|13ip . 



Q 


7T + 7T 

->7r°7r 


TT-R+ 

->ir°K° 


7T p 

— > 7r°n 


K+K- 


K-p 
-> K°n 


pp 

— » fin 


k*p, MeV/c 
fm 


35.5 

5.6 

-0.84 


11.3 
17.5 
-0.07 


28.0 

7.0 

-0.46 


i 62.9 
-i 3.1 
0.20 


i 58.6 
-i 3.4 
0.13 


i 49.3 
-i 4.0 
0.03 



k* — > — k* (the scattering is viewed in the diagram from right to left so that the final- 
state particles 1, 2 are in the entrance scattering channel). We will denote the channels 
as a = {1 + 2} and (3 = {3 + 4}, and the corresponding wave functions describing the 
scattering a — > a and a — > [3 - as ip a and ip 13 respectively. Outside the range of the 
strong FSI, r* > d, they can be written, for the a- and /3-channel continuous spectrum, 
as 



e- ik * r *F(-i^, 1, iU + 



V4*(r*) =A/^)r(F),/^%M (130) 



where A/"(r/ Q ) = e iSo ^ a ^ A c (r] 




[m 4 2 - m 3 2 + K + ^2*)T 

7/t4 

(131) 



2 _ L //6 4 — 1/1,3 -r yu^ -r U7 2 ; j 2 



= — fc* 2 + 2/i B (mi + m 2 -m 3 - m 4 ). 

The approximate equality in (j!3ip corresponds to the non-relativistic expressions for the 
energies: u* = mj+k^/ (2rrij),j = 1, 2. We consider here the systems with the Coulomb 
interaction absent in the channel j3, so ap = 00, r\p = 0, A c (r]p) = 1, G(pp, rjp) = exp(ipp) 
and x{Vp)/ a P = i^gj the amplitude ip 13 in ( 11301) then reduces to the expression indicated 
by the arrow. The /3-channel momenta at the a-channel thresholds (k* a = 0) for tttt-, ttK- 
, tcN-, KK-, KN- and iVTV-systems are given in table [3l This table also demonstrates 
that even close to the a-channel threshold, the use of the non-relativistic approximation 
can lead to noticeable shifts in kg. 



Similar to equation (1961) in the single-channel case, the amplitudes 

f xx ' = f xx '[A c (r ]x )M Vx/ )}- 1 / 2 , (132) 

where f xx ' are the amplitudes of the low-energy s-wave scattering due to the short- 
range interaction renormalized by the long-range Coulomb forces, A, A' = a, (3. The 
time-reversal invariance requires f xx ' = f x x . It is convenient to consider the amplitudes 
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fc X ' and / AA ' as the elements of the symmetric matrices f c and / related by the matrix 
equation 

f(k*) = [A c m l/2 fc(k*)[A c m 1/2 - (133) 

The single- channel expression ( |97|) for the amplitude f c can then be rewritten in a matrix 
form: 

f c (k*) = [k- 1 - , (134) 

where a, fj, x(v) an d A c (r)) are diagonal matrices in the (a, /3)-channel representation, 
for example, [A c (^)]^/ = A c (r)\)5\\>. The symmetric matrix K has to be real for the 
energies above both thresholds due to the two-channel s-wave unitarity condition [25] 

9f/ = / + «fc/, (135) 

where the diagonal matrix k\\i = k^5\\>. Usually, the AT-matrix is real also for negative 
kinetic energies (provided sufficiently fast vanishing of the short-range potential with 
the distance), and so it can be expanded in even powers of k* = fc*, similar to (IHTj) or 
(188]) with the parameters substituted by the corresponding matrices (e.g., f — > f ). 

Since, in the cases of practical interest, the particles (pions, kaons, nucleons) in 
the channels a and (3 are members of the corresponding isotopic multiplets, one can 
assume the parameter matrices diagonal in the representation of the total isospin |48j. 
The elements of the parameter or K (K~ l ) matrices in the channel representation are 
then given by the corresponding isospin projections. Particularly, for a = {tt + tt~}, 
(3 = {7r°7r }, one has: 

f aa _ 2 AO) , 1 A2) f af3 _ f l3a _ V% , AO) f (2)s 

Jo — qJo "t qJo Jo —Jo — ~^~\Jo Jo ) 

3 3 3 ( 136 ) 

f pi3 _ 1 AO) 2 (2 ) 

Analogous relations, with the substitutions (0) — > (1/2) and (2) — > (3/2), take place 

for the channels a = {k^p, tt~K + , tt + K~}, (3 = {ir°n, ir°K , tt°K }. For the channels 
a = {K + K~ , K~p,pp}, (3 = {K°K°, K°n,nn}, one has 

fo a = fS fi = ^/o (0) + /o (1) ) ff = ft = -\U^ - /o (1) ), (137) 

where the parameters and have now positive imaginary parts due to the effective 
inclusion of the additional channels opened at the energies of the elastic thresholds 
(fc* = 0) in the reactions KK — > tttt, nr], KN — > ir'sA, 7r'sS, A^A^ — > mesons. 

Note that the use of the isospin relations (I136P and (11371) implies that the violation 
of isotopic invariance is solely associated with the Coulomb factors A c (r]j) (strongly 
deviating from unity at fc* < 2tt / 1 | ) and the mass differences between the members of 
the same multiplets (fc* ^ k%). These relations however neglect the direct violation of 
isotopic invariance of order O(f /a) due to the renormalization effect of the Coulomb 
interaction on the scattering lengths, usually leading to the shifts on the level of several 
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percent. Within this uncertainty, one can also use (11361) or (11371) directly for the elements 
of the matrices K' 1 or K. 

The difference between the channel momenta can be neglected sufficiently far from 
the threshold. Then, one can apply (I136P or (I137P to the amplitudes fjj/ in the absence 
of the Coulomb interaction and, switch on this interaction in a similar way as in the 
single-channel case [53J: 

/>*) = [AM} 1/2 [Mk*) + ik - 1 [AM] 1/2 - (138) 

One may note that (11361) . (I137P correspond to the two-dimensional unitary 
transformation f = U^/qU, Un = U22 = cosy?, U12 = —U21 = sin (p. Since it 
applies also to the do-matrix, one immediately arrives at the same transformation of the 
complete amplitude / in the case of absent Coulomb interaction and k — k*l. 



6.2. Discrete spectrum in the a-channel 

One can repeat the same arguments as for the single-channel case, starting from the 
general solution in (11041) with the substitution f c — > f" a . For a discrete energy level 
E n = —n n 2 / (2/i), the amplitude f c has to have a pole or, equivalently, det / c _1 (m n ) = 0. 
Following [5QJ and introducing the matrix 

(i^) AA ' = (K- 1 )^' - iSxxSwk}, (139) 

one can rewrite this requirement in a form of equation f II 1 3 f) modified by the substitution 
K(iK n ) — * A aa (iK n ) and thus, fix the discrete energy levels similar to (11141) : 

K n = K c n {l + 2A aa K c n [l + A aa K c n [<t>{n) - 1] + 0(a- 2 )] } (140) 



a\2 



K aa - ifc* det K \k* a (K^ , 

where a = a a . Since K{iK n ) = K(0)[1 + TrO(/ ^o(™0~ 2 )] and k*^{iK n ) = k%(Q)[l + 
0((nakp)~ 2 )], one can safely make the substitutions K(iK n ) — > fo = K(0) and 
kp(iK n ) — > kp(0) and write, with the relative errors 0(a -2 ) less than a fraction of per 
mil, 

AA , _ K xx ' -ik*pdetK5 xx >5 Xa ^ / AA ' - det foS xy 5 Xa 



particularly, 



(k*K pa ) 2 (k* f^) 2 

jo ja« jv'aa P — foot fPP ^ fBJO / 

~ K K 1 + (k^r ~ /o /o 1 + 



SA C 



kUK^) 2 . kucr 



(143) 



In (I142p and (11431) . k% simply denotes k%(0) or k%(iK n ). It can be seen from table [3] 
that k^ 1 represents a scale which is intermediate between the Bohr radius \a\ and the 
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elements of the matrix f . As a result, the terms like 0(H(/q A ') 2 /a) or 0((ak*p)~ 2 ) 
contribute less than a fraction of per mil and can be omitted. As for the terms 
0((A;S/o A ) 2 )j their contribution is on a per mil level and is retained. 

The s-wave solutions corresponding to the a-channel discrete spectrum are again 
given by the second term in (110411 (A/"(r] n ) = 0) with the finite normalization N' = 
N fc° 7 ' A aa introduced in the same way as in fll 17j) modified by the substitution 
K — > A aa . As for the corresponding /3-channel s-wave solutions ipno( r *)> they are given 
by the second of equations ( 11301) with 

f/3a TS-0C 
.Jc ■ " J1 



ftffc = M'A C 



M'- 



M'A P ° 



(144) 



fr ' 1 

the second equality following from (1141jl and the explicit inversion of the symmetric 
matrix f~ lm . 



K c 



ik* det K 



det K = K aa K^ - [K^f 
\u\ 

D = det f- 1 det K = 1 - ihlKW + %{K aa - ikl det K), 



KM + ^ det if 



(145) 



where \ denotes here xiVa)] recall that x(.Vp)/ a i3 = ikt due to the absent Coulomb 
interaction in the channel f3. Note that the product Df xx> is finite since the amplitude 
pole for a bound state is compensated by the corresponding zero of the factor D oc 
det/" 1 . For the continuous spectrum at the a-channel threshold, K = f and 
D — 1 — iktC - (2i7r/|a|)(/ aa - ik% det /„). 



As a result, 

W) = M'{n) A< 



G 



n 



fr-ikldetf G 



1 - iAfc/ff 



1 + TrO 



o«o 



\n 2 a? / 




(146) 



1 + TrOl 



fpdp 
n 2 a? 



where G = G(p n ,r] n ) with the arguments p n and ^ n taken at fc* = i/« n (k„ is expressed 
through f in <[HD|> and (jI35jl ). and Af'(n) = ^^(O^l + OlfifV )] is fixed b Y the 
normalization integral ( 11 181) for the wave function ip^o- ^ can be calculated also 
analytically using (1 1 2 1 [1 with the substitution / — > ^,A aa rs /q"* (see Appendix B): 

2N 



|AA'(n)M°o ul (O)r-l = 0H 



no 



4tt 2 



UA C 



(147) 
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Using (11451) . one can express the amplitudes f xx '(k*) at k* a = (x = ~ m > K — fo) 
through the elements of the A-matrix (related to the scattering lengths /q A ' in (11421) ) 
with the relative error 0(a -2 ) less than a fraction of per mil: 

2m det f 



f^(0)=A 



XX' 



2l7T 

1 + --A aa + 0(a- 2 ) 

(X 



:i48) 



6.3. Universality 



Comparing (I130p and (I146p . one may see that the universal r*-behaviour of the s-wave 
amplitudes ^ A in continuous (k* —>■ 0) and discrete spectrum takes place with similar 
accuracy as in the single-channel case. Thus, using the expansions (I106p and (11191) 
for the amplitudes ip a , modified by the substitutions f — > /" a (0) and f — > A aa 
respectively, one has for the measures of the universality violation defined as in (11221) . 



+ 4vr"0 



^nO 



+TrO| 



TrO 



' fodp \ 
n 2 a? j 

( fod s 
{ n 2 a? 



+ 



+ 




+ 0(p 2 ) 



+ 0(p 2 



(149) 



The presence of the second channel manifests itself through a new scale hp (see table [3]), 
basically leading to the additional correction of 47rO(A;S(/ ( ^ a! ) 2 /a)) which is still on the 
negligible level less than a fraction of per mil. 

For the production cross sections, instead of (1751) and (T75|) . we now have: 

dV (V'a:: ^ . / ,.. . -'\ 



7i72 



d 3 pid 3 p 2 

+7374 



r\ 6 rr a 

u °0 \ ' n 

:7l72 dV^? s ' Q 

d 3 p 3 d 3 p 4 Y \ 



-k* ( 



qPS 



qPS 



(150) 



lb 



dVf 
d 3 P h 



(2tt) 3 7 1 7 2 



d 6 a£ 



-(27r) a 7 3 74 



d 3 pid 3 p 2 

d 6 ^ 



S,a, 



d 3 p 3 d 3 p 4 



OPS 



OPS 



(151) 



where p^ = P^m^/ (mi + m 2 ) in (11511) and b = {nO}. Since the particles 1, 3 and 2, 4 are 
usually the members of the same isospin multiplets, we can take 7i72d 6 o"Q = 7374d 6 o"o 
as a common factor in (11501) and (I15ip and also put Qs,a — Gs,p- 

The two-channel effects in the production cross section, being quadratic in the 
amplitude usually represent less than several percent of the strong FSI contribution 
(a fraction of percent in the cross section). Thus, for a near-threshold two-pion system 
produced according to a Gaussian redistribution ( 1159P with the characteristic radius 
rg = 3 fm and, taking the two-pion s-wave amplitudes from [8] (/o* a = 0.186 fm, 
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f$ a = —0.176 fm), the contributions of the FSI transitions 7r + 7r~ <-> 7r + 7r~ and 
7r + 7r~ <-> 7T°7r to the 7r + 7r~ production cross section respectively compose 7.72% and 
0.16%, ; these contributions are somewhat higher, 9.66% and 0.20%, for the amplitudes 
from [5] (f^ = 0.232 fm, Jq = —0.192 fm). At large rg, the elastic and inelastic 
contributions vanish as f$ a /rg and |/o / r </| 2 respectively. One should also account for 
the correction due to the deviation of the solutions in (I130p and (11461) from the exact 
ones in the inner region r* < d. Though this correction vanishes as r^ 3 , at rg = 3 fm 
it is still comparable to the contribution of the inelastic two-pion transition, composing 
0.25% and 0.20% for the amplitudes from [5J and [5j, respectively. 

Note that assuming 7i72d 6 <jQ Qs,a = Izlid^a^Qs^, the correction to the correlation 
function at a given total spin S, total four-momentum P and a small generalized relative 
four-momentum q = {0, 2k*} — > can be written as 



ATI 



d 3 r*WVr*) 



C k *(r*) 



+ 



+ 



(152) 



where Wp(r*) = j dt*gp(t*, r*; 0)/ J d A xgp(x; 0) is the normalized distribution of the 
vector r* of the relative distances between the emission points in the pair cm. system 
and ip denotes the solutions in (I130p extended to the inner region r* < d. In the case 
of only two open channels a and /?, the leading part of the correction scaled by Wp(0) 
is expressed through bilinear products of the amplitudes f xx ' in equation (44) of [48J. 
After a straightforward though lengthy algebra, it can be written in a more explicit 
form: 

caai2 d 



An 



+\rt 



47tW p (0)A c ( V , 
d 



-l\aa 



(3a\2_ 



dk 



*2 



dk* 2 



dk 



+2 



-l\/3a 



(153) 



at k* = 0, twice the derivatives of the inverse i^-matrix elements coincide with the 
effective radii d$ x ' . Similarly, in the case of discrete spectrum, the leading correction to 
(|^ (r*)| 2 + |V^ ( r *)| 2 ) is also S iven b y (HS3]) with the substitutions A c (r] a ) -> N'(n) 
and f xx ' — > A xx ' . For the Gaussian r*-distribution, (11531) is valid up to subleading 
contributions 0(k* 2 ai/rg) (see a discussion after (1581 ) and O^f^d^/rg). 

It is important that the presence of the second channel does not practically modify 
the ratio fl 1 2 5 [) of the finite-size correction factors in discrete and continuous (k* — > 0) 
spectrum at moderate distances r* <C \a\. The only modifications are the substitution 
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f — > 3?A Qa ps and the appearance of the negligible correction 47rO(A;^(/ ( 
l + 5 n _ [^(l^o(r*)| 2 )^ + ^(|^ (^)| 2 )^]|^r(0)r 2 



/3W0 ) 



! + *(**) [^(l^(r*)i 2 }g, + ^(i^(r*)| 2 }g J ]i« 1 ^ 1 - 2 




l + <J(ife*) 



1 + 



25M C 



n 



n a 



+ 0(F 2 (r* 2 ) SL ) +4ttO 



/M): 



(154) 



7. Finite-size effect in the experiment DIRAC 

7.1. 7r + 7r~ system 

We will use the results of the UrQMD transport code simulations of the pion production 
in pNi interactions at 24 GeV in the conditions of the DIRAC experiment at CERN 
[29] . Since we are interested in the region of very small relative momenta Q = 2k* < 20 
MeV/c, where the angular distribution of the vector Q is isotropic for non-correlated 
pions and, for Q < 10 MeV/c, the detector acceptance is practically independent of the 
direction of the vector Q, one can simplify the analysis integrating over this direction. 
The finite-size effect is then determined by the distribution of the relative distance r* 
between the pion production points in the pair cm. system, irrespective of the angular 
distribution of the vector r*. In fact, due to the applied cut Qt < Qt 11 = 4 MeV/c [I], 
this is true for Q < Qp* only, see the discussion after (199]) and (11011) ). For larger Q- 
values or, in the case of a two-dimensional (Qt, QL)-analysis, one needs two-dimensional 
(r T ,r1)- or (r*, cos ^^-distributions. We will neglect this complication here. 

The simulated r*-distribution is shown in figure [2J The tail of this distribution 
(r* > 50/m) is dominated by pion pairs containing a pion from the decays of to and 
rj' resonances, except for the pairs with both pions from one and the same decay. The 
respective decay lengths in the rest frame of the decay pion are about 30 and 900 fm; 
the decay length I ps r(pdec) /m^ is determined by the resonance lifetime r and the four- 
velocity Pdec/ m ir of the decay pion. As a consequence of the exponential decay law, the 
form of the corresponding redistributions is nearly exponential, except for the region 
of small r* dominated by the phase space suppression factor oc r* 2 . The exponential 
form is also distorted due to the averaging over the continuous spectrum of the decay 
momenta and over the emission points of the second pion. For r* less than 2000-3000 fm, 
the simulated r/ contribution (~ 1% of pion pairs at Q < 50 MeV/c) can be sufficiently 
well parametrized by an exponential-like formula interpolating between the phase space 
and exponential behaviour: 

2^ — ^ — = «y^( r I r rfM ( 155 ) 
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x 2 



F^lTrfJrf) = — |1 -exp 

& = — y = v 




1 + 0.2a; z 



1 + x 5 /125 



(156) 



where r^/ = 2 fm, = 790 fm. At the same time, a good description of the ui 
contribution (~ 19% of low-Q pion pairs) requires a superposition of two exponential-like 
expressions: 

V] = niuFir*; r^, l lw ) + n 2u T{r*\ r 2uJ , hu>)- (157) 

~r, ar* 

The parameters r Xw = 1.07 fm, = 43.0 fm, r 2ul = 2.65 fm, = 25.5 fm, 
niui/n2uj = 0.991 in the interval 2-200 fm and ri w = 1.00 fm, = 44.0 fm, r 2u) = 2.55 
fm, — 25.8 fm ni^/n^ = 0.845 in the interval 2-350 fm. We will use the former 
parameter set, but we have checked that the use of the latter one leads to a negligible 
change (< 0.1%) of the breakup probability. The rest of the redistribution due to 
the pions produced directly in the collision, in the rescatterings or in the decays of 
resonances with the decay lengths shorter than l u is peaked at ~ 3 fm and its main part 
(~ 60% of low-Q pion pairs) including the tail for r* = 10 — 100 fm can be effectively 
described by a power-like expression: 



M(r*;r M ,a,P) = r 



*2 



1 + 



2a 



-2/3 



;i58) 



r* V 

where r^ = 9.20 fm, a = 0.656, f3 = 2.86; note that the tail vanishes as (r*) -5 ' 5 , i.e. 
much faster than the Lorentzian (a = (3 = 1). The remaining short- distance part of 
the redistribution (~ 20% of \ow-Q pion pairs) is strongly shifted towards the origin 
because the UrQMD code assumes the point-like regions of the decays and rescatterings; 
particularly, r* = for ~ 8% of low-Q 7r + 7r~ pairs. Therefore, we will represent this 
part by a Gaussian distribution: 



g(r*;rg)=r* 2 exp[- , (159) 




where the Gaussian radius rg w 1 — 2 fm. As a result 



E ^T j) = n M M{r*- r M , a, 0) + n Q Q{r*- r g ). (160) 

We will also represent the short- distance part of the redistribution by the Gaussian 
contribution alone, i.e. put = and rg = 3 and 2 fm in (11601 ). 

The correction factors 1 + 5(k*) and 1 + 5 n corresponding to the redistributions 
r]',u),J\A,G, required to calculate the 7r + 7r~ production cross section in the continuous 
and discrete spectrum, are shown in figure [TlJ The two sets of histograms denoted 
by the same lines (dotted, full, dash-dotted, dashed and full) correspond to the two- 
pion scattering amplitudes from [8] (lower) and [5] (upper). In increasing order, 
they correspond to the redistributions r/, u, (?(r*;3fm), A4(r*; 9. 20fm, 0.656, 2.86) 
and (?(r*;2fm). One may see that the correction factors corresponding to the r\ 
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contribution are practically independent of the two-pion scattering amplitudes and, 
noticeably deviate from the infinite-size correction factors 1 + S°°(k*) = 1/A c (rj) (the 
curve) and 1 + 5^° = 0. We thus do not include the ?/-meson in the class of LL emitters, 
unlike the ?7-meson with the decay length of ~ 10 5 fm. 

The calculation of the correction factors was done according to the two-channel 
expressions given in the numerator and denominator of the first equality in fl 1 54j) : 



HL 



1 + 5(k*) = (|C k *(r*)f + |^ k *(r*) ~)„ p [Mv))-\ (161) 

/ 2\ SL —2 

l + 5n=(\W)\ 2 + \A(r*) \) k°o ul (0) , (162) 

\ /OP 

where a and (3 respectively denote the channels tt + tt~ and 7r°7r°. However, the account 
of the coupled tt 7t° channel and of the leading correction due to the approximate 
treatment of the wave function inside the range of the strong interaction, does not 
practically influence the results corresponding to the r] and to contributions and only 
slightly (< 1%) shifts up the correction factors corresponding to the short-distance M. 
and Q ones. A shift of the correction factors can arise also from the uncertainty in the 
s-wave elastic 7r + 7r~ scattering length f . The shift due to ~ 20% difference of the two- 
pion scattering amplitudes from [8] (fo = 0.186 fm) and [5] (/o = 0.232 fm) is ~ 2 — 3% 
for the short- distance M. and Q contributions and ~ 1% for the uj one. The global shifts 
are however not important since they can be absorbed in the product Xg = A in ( |33l ) 
and flS}. 

In accordance with the results in table [2] and figure [9] one may see in figure [11] the 
nearly universal slope of the factors 1 + S(k*) corresponding to the short-distance M. 
and Q contributions. In accordance with fllOip . the slope scales with f and is ~ 20% 
steeper when using the two-pion amplitudes from [5] instead of those from [8]. This is 
clearly seen in figure [121 where we plot the same correction factors as in figure [11] in 
a larger scale and with the subtracted intercepts 1 + 5(0). At Q > 20 MeV/c, there 
is also seen ~ 5 — 10% variation of the slope corresponding to different short-distance 
distributions. 

Figures [H] and [12] also demonstrate the violation of the universality relation 
5 n = 5(0) up to ~ 0.4% for the short-distance and uj contributions and up to ~ 9% for 
the rj' one. The most right panel in figure [12] shows that, in the case of the short-distance 
contribution, this violation is mainly related to the effect of the strong interaction on the 
normalization of the pionium wave function. Indeed, the difference 5' n — 5(0), corrected 
for this effect according to (EH), practically vanishes. 

In figure [13] we plot the correction factors corresponding to the mixture of 1% r/, 
19% uj and 80% short- distance contributions, as expected from the UrQMD simulation 
of low-Q pairs of charged pions in conditions of the experiment DIRAC. We neglect here 
the dependence of the contributions on Q, as well as the dependence on the pion charges. 
In fact, within the analysis region of Q < 15 MeV/c, the simulated u contribution 
increases with decreasing Q by ~ 0.01 and its average value for 7r + 7r~ pairs composes 
~ 0.15 [29]. 
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To show the effect of a possible uncertainty in the short-distance part, we describe 
it by the Gaussians with different characteristic radii rg = 3 and 2 fm. To account 
for the uncertainty in the two-pion scattering amplitudes, we have used those from [8] 
(/o — 0.186 fm) and [5] (/ = 0.232 fm). One may see that the corresponding global 
variations of the correction factors compose ~ 5% and ~ 2%, respectively. In figure [HI 
we plot the same factors with the subtracted values of the intercept 1 + 5(0). One may 
see that after the subtraction, the correction factors calculated for the same two-pion 
scattering amplitude but at different values of rg practically coincide for any n in discrete 
spectrum and for Q < 20 MeV/c in continuous spectrum. Since the subtraction can 
be included in the overall normalization factor, one may conclude that the uncertainty 
in the short-distance part of the r*-distribution is of minor importance for the relative 
momenta Q < 20 MeV/c. As for the effect of ~ 20% increase of the s-wave elastic 
7r + 7r~ scattering length, it leads to ~ 20% increase of 5 n — 5oo and to ~ 20% decrease of 
5(k*) - 5(0) at Q = 12 MeV/c. 

To estimate the effect of the uncertainties in the uj and rf contributions, we plot in 
figure[15]the differences 5—5(0), varying these contributions by ~ 30%. One may see that 
the corresponding variations of the differences respectively compose ~ 30% and ~ 20% 
for 5(k*) - 5(0) at Q = 12 MeV/c and, they are quite small (< 0.0003) for 5 n - 5(0). It 
should be noted that the rf contribution to the correction factor at Q > Q c is quite close 
to the infinite-size contribution 1/A C . The latter is included in the fit of the non-atomic 
7T + 7r _ correlation function thus essentially reducing the corresponding uncertainty in 
the breakup probability. This is demonstrated in table EJ where the contributions 
-AN A /N A and AN%/N% to the relative shifts AP br /P br = -AN A /N A + AN^/N^ 
of the breakup probability due to the neglect of finite-size corrections, corresponding 
to different mixtures of the rf, uj and Q contributions and different fit and signal 
intervals, are presented (see (I4"3"|) - (|4"T}) ). One may see that the 30% uncertainty in the rf 
contribution leads to negligible variations in the relative shifts AN A /N A (< 0.03%) and 
AN%/N% (< 0.2%). 

As for the uncertainty of the short-distance part of the redistribution, introduced 
by 30% decrease of the Gaussian radius from 3 to 2 fm, it also leads to negligible changes 
of the relative shifts AN A /N A (< 0.05%) and AN h A r /N h A r (< 0.4%) that rapidly decrease 
with decreasing upper boundaries of the fit and signal intervals. 

One can also neglect the present ~ 5% uncertainty in the 7r + 7r~ scattering length 
fo- Thus even the variation of fo by 20% leads to rather small variations of the relative 
shifts AN A /N A and AN^/N^ 1 ; e.g., for Q cut = 4 MeV/c and the fit interval (4, 15) 
MeV/c they compose only 0.12% and 0.96%, respectively. 

The dominant uncertainty in the finite-size correction to the breakup probability 
arises from the uncertainty in the u contribution. One may see from figure [16] that the 
correction APbr almost linearly increases with the fraction f u of 7r + 7r~ pairs containing 
a pion from u decay and the other pion from any short-lived source, except for pion 
pairs from one and the same uj decay; APbr/Pbr ~ —0.032 + 0.41/^ for Q C ut = 4 MeV/c 
and the fit interval (4, 15) MeV/c. Thus, taking f u = 0.19, a 30% (±0.06) variation in 
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Table 4. The contributions -AN A /N A and AN\ T /N\ Y to the relative shift 
AP br /P hT = -ANa/Na + AN^'/N^ of the breakup probability (U) due to the neglect 
of finite-size corrections. The non-atomic 7r + 7r~ correlation functions, calculated 
according to (|44j) for different mixtures of the rf t to and Q{rg) contributions, were 
fitted by ([45]) (fits i = 1 - 7) and AN A /N A and AN^/N^ 1 were calculated according 
to (f46|) and (|47|) . In approximate correspondence with [4], a uniform population of 
non-correlated pion pairs in Q was assumed in the considered fit intervals {Q\,Q2) 
and the ratios /N^ in the signal intervals (0, Q C ut) were set equal to 16, 9 and 

4 for Q cu t = 4, 3 and 2 MeV/c, respectively. 



FIT i 1 2 3 4 5 6 7 fit signal 

rg fm 3 3 3 3 3 3 2 region cut 

w% 19 25 13 19 19 19 19 Qi,Q 2 Qcut 

rf % 1.0 1.0 1.0 1.3 0.7 1.0 1.0 MeV/c MeV/c 

phase shifts [8] [8] [8] [8] [8] [5] [8] 

-AN a /Na % 1.09 1.61 0.58 1.11 1.06 0.89 1.14 4, 20 

AN^/N^ 1 % 4.06 6.77 1.30 4.22 3.86 2.69 4.48 4 

3.34 5.46 1.18 3.47 3.18 2.26 3.63 3 

2.08 3.37 0.80 2.17 1.98 1.44 2.24 2 



-AN a /Na % 1.02 1.43 0.61 1.04 0.99 0.90 1.03 4, 15 

AN^/N^ 1 % 3.62 5.57 1.50 3.70 3.42 2.66 3.82 4 

3.00 4.59 1.32 3.09 2.85 2.25 3.14 3 

1.89 2.88 0.88 1.96 1.80 1.44 1.96 2 



~AN A /N A % 0.83 1.11 0.55 0.85 0.81 0.78 0.83 4, 10 

AN^/N^' % 2.43 3.62 1.16 2.53 2.35 1.92 2.51 4 

2.12 3.14 1.06 2.21 2.04 1.70 2.17 3 

1.39 2.05 0.72 1.46 1.33 1.12 1.41 2 



f u leads to ~ 50% (±0.024) variation in AP br /P br . 

The correction rapidly decreases with decreasing upper boundaries of the fit and 
signal intervals. A decrease of the boundaries is however limited due to the increase of 
statistical errors. Also, the decrease of Q cut below 3 MeV/c introduces a systematic shift 
of ~ 5% in the breakup probability due to possibly insufficiently accurate description of 
the shape of the Q-spectrum of the atomic breakup 7r + 7r~ pairs [4J. The optimal choice 
seems to be Q cut = 4 MeV/c and the fit interval (4, 15) MeV/c. 

Taking f w = 0.15, the overestimation of the breakup probability in these fit and 
signal intervals composes 3% and corresponds to ~ 7.5% overestimation of the pionium 
lifetime. Correcting for this overestimation and assuming rather conservative 30% 
uncertainty in the to contribution, the uncertainty in the breakup probability composes 
2%, corresponding to ~ 5% uncertainty in the extracted pionium lifetime. 

One may expect an increase of the finite-size correction and its uncertainty when 
taking into account the increased slope of S(k*) at Q > = 4 MeV/c. On the other 
hand, one may expect a reduction of the correction when extending the fit region down to 
Q = and performing a more constrained 2-dimensional fit of the (Qt, QzJ-distribution 
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taking into account also the shape of the spectrum of atomic pion pairs. 
7.2. 7t~tt~ and 7r + 7r + systems 

As a by-product, the experiment DIRAC provides data on the correlation functions of 
identical charged pions which contains the information on the space-time characteristics 
of pion production and can be used to check the results of the UrQMD simulations. 
Thus, ~ 5 ■ 10 5 7r~7r~ pairs have been collected in pNi interactions in 2001, representing 
~ 40% of the available statistics [29J. The Q-distribution of these pairs is peaked at 
~ 60 MeV/c and drops essentially outside the interval (20, 120) MeV/c due to a decrease 
of the phase space and detector acceptance at small and large Q respectively. 

Contrary to the case of the 7r + 7r~ system, the correlation effect in the system of 
identical pions extends and is measured up to the relative momenta Q ~ 200 MeV/c, so 
neither the distribution of the vector Q nor the detector acceptance can be considered 
independent of the direction of this vector. Since further the angular distribution of the 
vector r* is not isotropic (particularly, the characteristic width of the out component 
of the redistribution increases with the transverse momentum while those of the side 
and longitudinal ones decrease), the required space-time information does not reduce to 
the distribution of the relative distance r* between the pion production points in the 
pair cm. system; generally, the 3-dimensional distribution of the vector r* is required. 
Here we however neglect this complication and calculate the 1-dimensional correlation 
function of two identical charged pions in the same way as for the previously considered 
case of the near-threshold 7r + 7r~ system, i.e. assuming the uniform distribution of the 
cosine of the angle between the vectors Q and r* for the non-correlated pions. 

The calculated 7r~7r~ correlation functions 7l~r , TZ~~, H~m and 7Zg~ corresponding 
to the redistributions rj', u, M(r*; 9.20fm, 0.656, 2.86), and G(r*;r g ), r g = 3,2,1.5 fm, 
are shown in figure [T71 In figure [181 we show the correlation function corresponding 
to 1% r}', 19% u, 60% A4(r*;9.20fm, 0.656, 2.86) and 20% Q(r*; 1.5fm) contributions, 
as expected for the pairs of charged pions from the UrQMD simulations; the errors are 
taken from the DIRAC pNi 2001 data. To demonstrate the sensitivity to the relative 
uj contribution, / w , we show in this figure also the correlation functions calculated with 
f u varied by ~ 30%. One can conclude, that the different shape of the u contribution 
as compared with the shapes of the short-distance ones (A4 and Q) allows, in principle, 
to determine f w - the most critical parameter required to calculate the finite-size 7r + 7r~ 
correction factors. To estimate the statistics required to determine f u better than to 
30%, we have fitted the correlation function in figure [TBI by 



The parameter N cares for possible normalization mismatch, the correlation strength 
parameter A takes into account the contribution of LL emitters, particle misidentification 
as well as coherence effects and the slope parameter b cares for a possible mismatch in 



K{Q) = N{\ [UKj-(Q)+fJl—(Q)+fgKj-(Q;rg) 
+ (1 - U -f»- fg)n^{Q- r M , a, (3)] + (1 - A)} (1 + bQ). 



(163) 
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Table 5. The results of the fits according to (|163p of the tt~it~ correlation function 
corresponding to the mixture of the 1% 77', 19% ui, 60% M(rj^,a, (3) and 20% Q{rg) 
contributions with rg = 1.5 fm, r» = 9.2 fm, a = 0.656, (3 — 2.863. The errors 
correspond to the statistics of 5 TO 7 7r _ 7r _ pairs from pNi interactions in the conditions 
of the experiment DIRAC. Only the parabolic errors are shown, absent error means a 



fixed parameter. 


FIT i 


1 


2 


3 


4 


5 


6 


N 


1.000 


0.999 


1.000 


1.002 


0.998 


1.017 




±0.001 


±0.001 


±0.001 


±0.001 


±0.001 


±0.001 


A 


1.001 


0.993 


0.999 


0.973 


1.024 


0.984 




±0.004 


±0.003 


±0.004 


±0.003 


±0.003 


±0.003 


r M fm 


9.253 


9.239 


9.161 


7.381 


10.844 


6.322 




±0.048 


±0.049 


±0.048 


±0.045 


±0.065 


±0.032 


a 


0.657 


0.651 


0.655 


0.611 


0.746 


0.561 




±0.003 


±0.003 


±0.003 


±0.002 


±0.004 


±0.001 




2.864 


2.893 


2.853 


2.717 


2.813 


3.001 




±0.010 


±0.010 


±0.009 


±0.008 


±0.013 


±0.006 


rg fm 


1.495 


1.492 


1.496 


1.505 


1.554 


1.5 




±0.016 


±0.016 


±0.015 


±0.022 


±0.013 




fs 


0.202 


0.199 


0.200 


0.147 


0.261 


0. 




±0.003 


±0.003 


±0.003 


±0.004 


±0.003 




U 


0.187 


0.190 


0.184 


0.131 


0.238 


0.234 




±0.004 


±0.005 


±0.005 






±0.004 


u> 


0.010 


0. 


0.01 


0.01 


0.01 


0.01 




±0.004 












b (GeV/c) 


- 1 0.001 


0.000 


0.000 


-0.033 


0.036 


-0.162 




±0.008 


±0.008 


±0.007 


±0.008 


±0.008 


±0.006 


x 2 


0.01 


0.02 


0.01 


1.01 


1.01 


33.9 



the Q- dependence of the reference sample in the denominator of the correlation function. 
Equation (11631) does not take into account a decrease of the fractions f v ' and f u with 
increasing Q which can amount to ~ 20% in the interval 0-200 MeV/c [29|. 

The dependence of the correlation function on the parameters rg, r^, (3 is 
calculated with the help of quadratic interpolation. For example, to calculate the rg- 
dependence, three values of this parameter are chosen and the corresponding correlation 
functions TZg'iQ] r l g ), i = 1, 2, 3, are used to interpolate to the correlation function at a 
given value of rg according to (!82|) with the substitutions 7?.rqmd - ► Tig an d s r — > rg. 
Similarly, the quadratic interpolation in the three-parameter space rx, a, (3 requires 
the calculation of 3 3 correlation functions 1ZJ^~(Q; r l M , a\ /3 k ), i,j, k = 1, 2, 3. 

The fit recovers the input parameters with rather small parabolic errors, 
particularly, / w = 0.189 ±0.048. It appears, however, that the real error, corresponding 
to the increase of the x 2 by one unit, is one order of magnitude larger. Thus, to achieve 
the error in f u smaller than 30%, one has to collect the statistics of ~ 5 • 10 7 tt~tt~ 
pairs. The fit results corresponding to such a statistics are shown in table El fits 1-6. 
Comparing fits 1-3, one may see that they are practically insensitive to the nearly flat 
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rf contribution to the correlation function, except for ~ 20% drop in the first 5 MeV/c 
bin. Since the drop contribution of 0.2 ■ 0.01 is much smaller than the correlation 
function error of 0.02 in the first bin, the 0.01 shift in /,,/ in fit 2 (f v i = 0) leaves the 
parameters and x 2 practically unchanged, except for the compensating ~ 0.01 shift in 
the correlation strength parameter A. Comparing further fits 3-5, one can conclude 
that the true error in f u , corresponding to the increase of x 2 by one unit, composes 
~ 25% (~ 0.05). Finally, fit 6 (fg = 0) shows that the oversimplified description of the 
short- distance contribution could lead to ~ 25% systematic shift of f w . 

To infer the fraction /+ _ for 7r + 7r~ pairs at Q — > from the fitted for identical 
charged pions, one has to take into account its (^"dependence as well as the fact that, 
due to a lower multiplicity of the pairs of identical charged pions, the fraction is 
~ 40% higher than f+~ [29]. 

8. Conclusions 

We have developed a practical formalism allowing one to quantify the effect of a 
finite space-time extent of particle emission region on the two-particle production in 
continuous and discrete spectrum. We have shown that one can usually neglect the 
non-equal emission times in the pair cm. system, the space-time coherence and the 
residual charge. The developed formalism is in the basis of the femtoscopy techniques 
allowing one to measure space-time characteristics of particle production as well as 
the low-energy strong interaction between specific particles. We have applied it to the 
problem of lifetime measurement of hadronic atoms produced by a high-energy beam 
in a thin target, particularly, to the measurement of pionium lifetime in the experiment 
DIRAC at CERN. Based on the transport code simulations, we have calculated so 
called correction factors that can be used to take into account the finite size of the 
production region by multiplying the point-like Coulomb production cross sections of 
the free and bound 7r + 7r~ pairs. We have shown that the short-distance contribution 
is of minor importance for the lifetime measurement since it leads to practically the 
same and nearly constant correction factors for free and bound pairs which cancel in the 
breakup probability. The most important is the fraction / w of tt + tt~ pairs containing a 
pion from uj decay and the other pion from any short-lived source, except for pion pairs 
from one and the same u decay. Besides leading to slightly different global shifts of the 
correction factors, it also affects their Q- and n-dependence. The resulting correction to 
the "point-like" pionium lifetime composes ~ —7.5%. Assuming rather conservative 30% 
uncertainty in / w , due to the uncertainty in the u production and a simplified treatment 
of the correction (e.g., the neglect of f u variation in the analyzed Q- interval), one arrives 
at ~ 5% uncertainty in the extracted pionium lifetime. It is shown that this uncertainty 
could be diminished if the high statistics data on correlations of identical charged pions 
were collected in the DIRAC experiment. The statistics required to determine f u to 
10% and control the finite-size effect on the lifetime to ~ 2% composes ~ 3 ■ 10 8 tt~tt~ 
pairs. The uncertainty in / w can be also reduced by tuning the transport simulations 
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with the help of experimental data on particle (resonance) spectra and femtoscopic 
correlations in proton-nucleus collisions at the beam energy of ~ 20 GeV. The lifetime 
uncertainty can be essentially reduced in future experiments using the multi-layer targets 
[2] since it will be basically determined by the uncertainty in the calculated number of 
produced atoms only; even for the conservative 30% uncertainty in f u , the corresponding 
uncertainty in the lifetime will be ~ 1% only. The above estimates of the finite-size 
correction to the pionium lifetime and its uncertainty are based on a 1-dimensional fit of 
the Q-distribution in the interval (4, 15) MeV/c. One may expect their underestimation 
due to the neglected increase of the slope of Q-distribution for Q > Qr* = ^ MeV/c. 
On the other hand, one may expect their reduction when extending the fit interval 
down to Q — and performing a more constrained 2-dimensional fit of the (Qt, Ql)~ 
distribution. The effect of the finite-size uncertainty on the breakup probability of other 
hadronic atoms remains to be studied. For ttK and up atoms, it is expected similar 
to the one for tt + it~ atoms since ~ 50% decrease of the u contribution is compensated 
by about the same decrease of the Bohr radius \a\ thus retaining a similar finite-size 
Coulomb FSI effect ~ (r*) SL /a. As for K + K~, K~p and pp atoms, there is no u 
contribution to the corresponding hadron pairs though, due to the smaller Bohr radii, 
its effect can be partly substituted by the contribution of 0-meson and other sufficiently 
narrow resonances. 
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Appendix A. Non-equal emission times 

We consider here the role of non-equal emission times in the Bethe-Salpeter amplitude 
i^( x ) = e 1 ^/ 2 + Aip~~\x), where the correction Aip to the plane wave is given in 
(flOl) . We will consider the amplitude in the pair cm. system, in which the plane wave 
e iqx/2 _ e -ik*r* j g j nc l e p enc [ en t of the emission times. First, we will prove the integral 
relation between the Bethe-Salpeter amplitude and the corresponding non-relativistic 
wave function, derived on the condition k* 2 <C /i 2 [12J: 

4>i + \x) = J dV<5 t .(r- - r',i>_ t .(r') (A.l) 

M"'-^*) - (^/d 3 --^-'exp(-i^| f l) , (A.2) 

where m(t* > 0) = and m(t* < 0) = mi. 

We start by splitting the product of the propagators into four terms, each containing 
only two poles in the complex K -plane, situated in the opposite upper and lower half- 
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planes. Taking into account that in the pair cm. system P = and that the pair energy 
coincides with its effective mass: Pq = m 12 , we get 

{(k 2 - mi 2 + iO)[(P - k) 2 - m 2 2 + iO]}" 1 

= [(/c Q + iO)(/c + £1 - iO) 

•(/C - ?71i2 - ^2 + lO)(/C - 7Tii2 + u) 2 - iO)]" 1 

= [m 12 2 - - c^TMt^o - Si + iO)(/c + a>i - iO)]- 1 (A.3) 

+ [(«0 - "7,12 - ^2 + i0)(« - "7.12 + ^2 ~ i0)] _1 

—[(/Co — a)i + i0)(«o - 7?7i2 + ^2 - iO)] -1 
-[(/c + ui - iO)(/c - 777i2 - w 2 + iO)] -1 }, 

where u)i = {m 2 + k 2 ) 1 / 2 . Assuming now that the amplitude f s = / s (/co,77ii2 — /Co) 
is an analytical function in the complex /co-plane, we can integrate over k q using the 
residue theorem. Consider first t* > 0. In this case the integration contour has to be 
closed in the upper half-plane, equation ffTOl) then giving 

1 _;r™ l/-™. 2 ,_2-w, u* /o /" d 3 Ke~ lK,r 



A^l +) (x) = -^m 12 e- i[mi2+(mi2 - m22)/mi2r/2 / 

9 7T 2 J 



77712 2 - (^1 - U2) 2 



-iujlt* 



1 1 



e -/(-^mia + d) - (A.4) 

V 77712 + Wi + W 2 ^1/ 

Vr77i 2 — co>i — cj 2 + iO 2co>2/_ 
Since we are interested in the limit of small particle momenta in the pair cm. system: 
k* 2 <C fi 2 and since the integral (1A.4j) is dominated by k, 2 « k* 2 , we can use the following 
non-relativistic approximations (recall that /J = mxm^j (mi + 777 2 ) is the reduced mass 
of the two-particle system): 

k* 2 

77712 = 777l + 777 2 + — - 

2/i 

9 2\ / / 777 2 k* 2 \ 

7T712 + (7771 ~ 777 2 )/777l2 = 2 OTl H ; — 

V 777i + 777 2 2/i / 



K 2 

= 777j + m i2 2 - (£1 - = 4mim 2 



(A.5) 



c* 2 - k 2 



77712 - UJ X — UJ 2 



2/i 

Retaining in the integral (1A.4I) only the dominant pole term ~ [77712 — uj\ — uj 2 + iO] -1 , 
we get 

A ,(+), s 1 f d 3 Ke- iKr * 



9 v ' 2tt 2 7 k 2 - k* 2 - iO 



2 ,,* 2 x (A-6) 



• expf-i K 2 ^ k t* j /(77712 - 2 2 , ^2 
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(+), 



J d 3 At'5 (3) (K- k' 



d 3 



ree 




e -iK(r*-r') exp 




(A.7) 



f{m 12 



-i • 



CJ 2 , L^o 



2tt 2 J K /2 -k* 2 -i0" 
iV4»(r*-r',r)A^ k ,(r'), 
where the ^.-function is given in (IA.2I) . Noting that the ^.-function in the integral 



(1A.7j) acts on the plane wave e _lk * r ' as a 5-function, we finally arrive at the integral 
relation in (1A.1I) for t* > 0. The prove of this relation in the case of t* < is done 
in a similar way, the integration Ko-contour being now closed in the lower half-plane. 
The result is the same as in (1A.4I) - (1A.7I) . up to the substitutions m% — > —m,\ in the 
time-dependent phase factor and 0)2 — > raw — U)\ in the arguments of the scattering 
amplitude /. 



At t* 



0, the function 5k* (r* — r', 0) : 

/ m 2 \ 3 / 2 



5 k *(r*-r',t* 



V27rit =t 



5 {3) (r* -r') and, at t* > 0, 
(k*H* 



cxp 



+ 



r') 2 m 2 N 



(Ai 



2m 2 2t* 

For negative t*-values, the substitution m 2 — > —mi has to be done in flA.8[) . It is clear 
from (1A.8j) that, at small k* (k* m(t*)r* /\t*\), the function 5fc*(r* — r',t*) practically 
coincides with the 5- function <5®(r* — r') on condition (!73|) . 

Since the particles start to feel each other only after both of them are created, it 
is clear that a large difference in the emission times generally leads to a suppression of 



particle interaction at small k*\ \ Aip~~' (x) \ < |A^_k*(r*)|; Aij;~(x) — at /- : — ' x.. 



Particularly instructive is the case when one of the two particles is very heavy, say 
m 2 ^> mi. Then the two-particle interaction is suppressed provided the light particle 
is emitted prior the emission of the heavy one (m(t* < 0) = m\ in (IA.2I) ). Otherwise, 
the large mass m(t* > 0) = m 2 prevents the suppression even if the light particle were 
emitted much later than the heavy one. Below we consider the effect of non-equal 
emission times on two-particle production in some detail. 

We start with the FSI due to the short-range forces only. Inserting the spherical 



wave fl85l into the integral relation flA.ll) or (IA.7I) . we get [12 

A^l +) 



/(**) 



where z ± = (^) 1/2 (r*±™) and E 1 (z) 



i sin(A;*r*) + 



Ei(z_ 



+ Ei(z+)e 



-ik*r* 



(A.9) 



iJdye 



is the Fresnel integral. 



Note that the length k*\t*\/m(t*) = ly.* can be interpreted classically, for large k*r*, 
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as a distance traveled by the first emitted particle until the creation moment of the 
second one. The absolute value of the factor (r* ± 4*) in the argument z + thus 
corresponds to the maximal (minimal) possible distance between the particles in their 
cm. system at the later of the two creation moments. The effect of non-equal emission 
times however doesn't reduce to the modification of the distance r*, it survives even at 
k* = 0. This effect vanishes in the limit of small \t*\, when z_ ^> 1, Ei(z±) — > (1 + i)/2 
and (IA.9I) reduces to the spherical wave fl85l) . In the opposite limit of large \t*\, when 
[m(t*)r* 2 / (2\t*\)} 1 ^ 2 <C 1, the interaction is suppressed and the scattered wave Aip~~\x) 
tends to zero for arbitrary fc*-values. 

In the simple static Gaussian model of independent one-particle emitters described 
by the amplitude (15T1) . the applicability condition (1731) of the equal-time approximation 
can be roughly written in the form ( |74l) . Clearly, the latter condition is not satisfied for 
very slow particles emitted by the emitters of a long lifetime. This is demonstrated in 
figures [3] and H] for the FSI contribution in the 7r°7r° correlation function. 

Note that the change of the character of the effect of non-equal times at v ~ 0.6 
and its increase with the increasing velocity is not expected from condition (jTlj) . The 
increase of the effect for relativistic particles (v — > 1) is specific for the systems of not 
very large sizes and lifetimes r ~ r , when the population of the light-cone region 
r ~ vt is not negligible. Indeed, in this region the arguments of the Fresnel integrals at 
k* = can be small even at large 7: z± « (jm\ri — 1|/2) 1//2 , leading to the modification 
of the spherical wave. 

Consider finally the effect of non-equal emission times on the correlations of two 
charged particles. Since, at not very large \t*\, the function 5k* {r* — r is close to the 
5-function, we can neglect the terms of higher powers of (r'/a) in (lA.ip fjj The non-equal 
time correction is thus mainly generated by the subleading term ~ r*/a and so can be 
expected rather small, similar to the case of strong FSI, where it arises from a small 
finite-size contribution ~ //r*. It concerns also the case of hadronic atoms since the 
Schrodinger equation at a small negative energy — e& = — K 2 /(2/i) practically coincides 
with that in continuous spectrum at zero energy. As a result, for r* C k" 1 = n\a\ (n 
being the main atomic quantum number), the r*-dependence of the corresponding wave 
functions at a given orbital angular momentum is the same. 

Appendix B. Decay rate and normalization 

The decay rate (partial width) of a bound a-channel state decay into the /3-channel is 
given by the square of the wave function ip^o i* 1 flHBI) ( a t a distance r* > d), multiplied 
by the product of the surface 4irr* 2 and the relative velocity v@ = k%/ {J,p: 




(B.l) 



1 + (klKPf 3 ) 2 ' 



J The account of these terms is however important at large time separations to guarantee vanishing 
of the Coulomb interaction at \t*\ — > 00. 
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Note that for two identical bosons in the channel /3, the twice as large square of the 
symmetrized wave function is compensated by twice as small surface so that the result 
is the same as for two non-identical particles. In the considered two-channel case, the 
/3-channel is the only open one, so the decay rate coincides with the inverse lifetime 
(total width) of the bound a-channel state which can be calculated from the imaginary 
part of the energy E n = — K n 2 /(2/j la ): 

l/r n = T n = -2%E n ee 2^K n QK n //i a . (B.2) 

Using ( 11 40 p . one has (neglecting QA aa as compared with the ?R.A aa in the correction 
terms) 



n 

2QA aa 



1 + 29fL4 ac X + 0((2gfL4 Q X 



;<) 2 {l + 2[<j>(n) - l]9fL4»X - 4vr 2 0((^/a) 2 ) } 



(B-3) 



T n = — {l + 2(j){n)^A aa K c n - 47r 2 0((5LT7a) 2 ) } «) 3 $M ac \ (B.4) 



Using the relation (k^) 3 
with [501 [M]: 

k 



7t|^q u1 | 2 and (j!43p for ^sA aa , one finally gets, in agreement 



1 + 20(n) 




^A c 



n\a\ 



a\2 



(B.5) 



\ + {k}K^) 2 ' 

Inserting (IB. 1|) and (IB. 5|) into the equality T n = one proves the relation (11471) 
between the normalization factors M'{n) and "0^ (0). 

In the case of two or more open decay channels, the two-channel (a,/3) matrix K 
is no more real, particularly, in the presence of one additional channel j, one has: 

ik*{JCi a ) 2 



K c 



r + 



1 - itf&i 



K f3o 



1 - iA££" 



(B.6) 



1 - ikj&j 



where /C AA ' are the elements of a real three-channel matrix /C. Note that in the case 
of a two-pion system (a = 7r + 7r~, (5 = 7r°7r°), the third channel is j = 77 so that the 
elements K)^ and /C" can be safely neglected. Then, only the element K aa acquires the 
imaginary part: K aa = IC aa + ik*{K j(X ) 2 . 

Generally, one has to account for the possible imaginary parts of the elements of 
the two-channel .fT-matrix as well as, for a possibility of a pure imaginary value of the 
momentum kp in the case of a closed channel [3 {kf < 0, k* p = i(-kf) 1 ^ 2 , = 0). 
Then 

MA aa = §lK aa - {Qk^iK 130 ) 2 + m^K^) 2 



+ 



k 



a\2 



1 - \k*pKM 



(B.7) 
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+ 



h* 



K fi 



1 - ik* p KW 



1 - ikpKM 



J a 



1 - ikpKW 



(B.8) 



4vr|iV'(n)| 2 



n 1 / j n 



where 9{x) = 1 for x > 0, 



for x < 0. The second expression for QA 010 in ( IB. 81 ) 



follows from a straightforward though lengthy matrix algebra and, the last one - from 
an obvious generalization of ( IB. II) using the relation K^ a = Df^ a . Inserting the last 
equality in (jR8|) into (IB~4l) . one proves (I147P for the case of any number of open decay 
channels. 
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Figure 1. The diagrams describing production of particles 1 and 2 in continuous and 
discrete spectrum. 
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Figure 2. The distribution of the relative distance r* between the pion production 
points in the pair cm. system simulated with the UrQMD transport code [28] for pNi 
interactions at 24 GeV and the relative momenta in the pair cm. system Q = 2k* < 50 
MeV/c in the conditions of the DIRAC experiment at CERN [29]. The curves are the 
results of the fits to short-distance, oj and rf contributions described in the text. 
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Figure 3. The FSI contribution to the ir°ir° correlation function calculated for 
different values of the pair velocity v in a model of independent one-particle emitters 
distributed according to a Gaussian law with the spatial and time width parameters 
ro = 2 fm and To = 2 fm/c. The exact results (solid curves) are compared with those 
obtained in the equal-time approximation (dash curves). 
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Figure 4. The same as in figure [3] for the pair velocity v = 0.1, the spatial width 
parameter ro = 2 fm and different values of the time width parameter tq . 
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Figure 5. The ratios of the ir + ir~ correlation functions 1Z Z and 1Z" Z " . For the latter, 
only one-particle spectra are influenced by the effective comoving charge Z. The pions 
are assumed to be emitted in the rest frame of a point-like charge Z according to 
the thermal law with a temperature of 140 McV. The distribution of the space-time 
coordinates of the particle emitters is simulated as a product of Gauss functions with 
the equal spatial and time width parameters r — ctq. The full broken line corresponds 
to Z = 30, r = 2 fm, the dash and dotted ones - to Z — 60, r = 2 and 3 fm, 
respectively. 
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Figure 6. The functions A c (rf), ^xiv) = ^(v) an d ^x(v) = A c (r])/(2r]) defined in 
([24]) . ([94]) and (|95|) . The solid and dash curves correspond to the attraction (77 < 0) 
and repulsion (rj > 0) respectively. For two-pion systems, the variable l^l -1 = \ak*\ 
approximately coincides with the relative three-momentum Q = 2k* in MeV/c: 
l^l -1 = 0.98Q/(MeV/c). The arrow in the first panel indicates the characteristic 
width = 2tt of the Coulomb effect. 
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Figure 7. The functions B(p,rj) and P(p,rf) defined in (|92|) . (f93|) and calculated for 
the 7r + 7r - system. The solid, dash and dash-dotted curves correspond to r* = 5, 15 
and 50 fm respectively. The dotted curves represent the functions B(p,0) = sin p/p 
and P(p, 0) = cos p corresponding to the case of neutral particles. 
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Figure 8. The function G(p, rj) defined in (|9Tj) and calculated for the ir + ir system. 
The solid, dash and dash-dotted curves correspond to r* = 5, 15 and 50 fm respectively. 
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Figure 9. The 7r + 7r~ correlation function at a fixed separation r* divided by the 
Coulomb penetration factor: 1Z/A C = (|c~ lk r F + f c G/r*\ 2 ), and the corresponding 
main contributions due to the interference term and the modulus squared of the 
hypergeometric function (see (|99|) ). The solid, dash and dash-dotted curves correspond 
to r* — 5, 15 and 50 fm respectively. The calculation is done in the approximation 
of a constant scattering amplitude f c {k*) = fo = 0.232 fm, the averaging assumes the 
uniform distribution of the cosine of the angle between the vectors r* and k* = Q/2. 
The dotted curves in the lower panel represent the s-wave Coulomb contribution 
B 2 {p 1 rf) to the quadratic term. 
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Figure 10. Comparison of the approximate ir + ir~ atomic wave function V'no P ( r *) anc ^ 
the 7r + 7r~ wave function in continuous spectrum r/)_k*(r*) at k* — > 0, respectively 
defined in (|120[) and ([89]) (/o = 0.232 fm), with the exact s-wave solution 
outside the range of the strong interaction Vno( r *) given in (|116|) : AfP p (r*) = 
[lCV)/tMr*)] a -landA*"^(r*)H^ 

fe* — > 0; the averaging in the latter expression is done over the uniform distribution of 
the cosine of the angle between the vectors r* and k* = Q/2. The central panel shows 
A^ pp (r*) assuming TV' (n) = tp^iO) in (|120|) in correspondence with the ansatz (|125j) . 
The curves in the increasing order correspond to n = 1, 2, 3, 10. 
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Figure 11. The correction factors 1 + S(k*) (left) and 1 + S n (right) as functions of 
the relative momentum Q = 2k* and the main atomic quantum number n respectively. 
They are required to calculate the tt + it~ production cross sections in the continuous 
and discrete spectrum according to (|33|) and ([34"]) . The two sets of histograms denoted 
by the same lines (dotted, full, dash-dotted, dashed and full) correspond to the two- 
pion scattering amplitudes from ,8i (lower) and [5] (upper). In increasing order, they 
correspond to the r* -distributions n' , uj, £(?-*; 3fm), M(r*; 9. 20fm, 0.656, 2.86) and 
C/(r*;2fm) defined in (fT55jl - l(T59"ll . The calculation was done according to the two- 
channel expressions given in ()161|) and (I162j) . taking into account the correction in 
(|153p . Note that the infinite-size correction factors 1 + 5°°(k*) = 1/A c (rf) (the curve) 
and 1 + 52° = 0. 
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Figure 12. The differences S(k*) - (5(0) (left panel), S n - (5(0) (middle panel) 
and (5^ — 5(0) (right panel) calculated from the tt + tt~ correction factors given in 
figure [TTJ The latter difference is corrected for the effect of the strong interaction 
on the normalization of the pionium wave function according to (|4ip . The differences 
corresponding to the rf contribution (dotted histograms) are not seen except for the 
first bin in the left panel; in the middle and right panels they compose ~ —0.1. 
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Figure 13. The ir + ir~ correction factors l + 6(k*) (left) and l + 5 n (right) calculated 
in the same way as in figure [TT1 assuming the mixtures of 1% rj' , 19% lu and 80% Q 
contributions with rg — 3 fm (lower and middle) and rg = 2 fm (upper). The lower 
and upper histograms correspond to the two-pion scattering amplitudes from [8] , the 
middle one - to those from [5]. 
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Figure 14. The differences S(k*) — 6(0) (left) and 6 n — 6(0) (right) corresponding to 
the correction factors in figure [T3l 
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Figure 15. The differences 6{k*)-S(0) (left) and S n -S(0) (right). The full histogram 
coincides with that in figure [T4l The dashed ones correspond to the 0.19 ± 0.06 Li 
contributions and the dotted ones - to the 0.010 ± 0.003 rf contributions. 
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Figure 16. Relative shift of the breakup probability AP\, T /P\, r = —ANa/Na + 
AN^ 1 /Njf (full curve) due to the neglect of the finite-size effect as a function of the 
fraction f u of tt + tt~ pairs containing a pion from ui decay and the other pion from any 
short-lived source, except for pion pairs from one and the same uj decay. The fit and 
signal intervals are (4, 15) and (0,4) MeV/c respectively (see table [4]). Also shown are 
contributions of the relative shifts —ANa/Na (dotted curve) and AN^/N^ (dashed 
curve) in the calculated numbers of produced atoms and breakup atoms respectively. 
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Figure 17. The ir w correlation functions. The histograms in the increasing order 
of the peak values correspond to the redistributions rf , u), M(r*; 9.20fm, 0.656, 2.86), 
C?(r*;3fm), g(r*;2fm) and 0(r*;1.5fm) respectively. 
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Figure 18. The ir~Tr~ correlation functions. The middle histogram and the fit curve 
correspond to 1% rf, 19% u, 60% M(r*; 9.20fm, 0.656, 2.86) and 20% £(r*;1.5fm) 
contributions. The errors are taken from the D1RAC pNi 2001 data [29]. The upper 
and lower histograms correspond to the 1% rf , 19% ± 6% lu contributions and the 
unchanged ratio 3:1 of the M. and Q contributions (unchanged form of the short- 
distance contribution). 



